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1.  INTRODUCTION 


Recent  years  have  seen  a renewed  interest  in  urban 
rail  transit  systems  for  public  transportation.  In  part  this 
interest  has  come  about  because  of  the  many  environmental  prob- 
lems associated  with  other  non-rail  transit  systems.  However, 
environmental  problems  are  also  associated  with  rail  transit 
systems.  This  report  focuses  on  one  problem  of  particular  con- 
cern; vibration  and  low  frequency  noise  in  buildings  near 
underground  transit  lines.  A second  problem  of  concern,  noise 
near  elevated  rail  transit  structures,  is  dealt  with  in  another 
report  prepared  by  Cambridge  Collaborative,  Inc.,  for  the 
Department  of  Transportation  [1] . 

1 . 1 Background 

As  systems  manager  for  the  U.S.  Urban  Mass  Transporta- 
tion Administration  (UMTA)  Rail  Systems  Supporting  Technology 
Program,  the  Transportation  Systems  Center  (TSC)  is  conducting 
research,  development  and  demonstration  efforts  directed  towards 
the  introduction  of  improved  technology  in  urban  rail  systems 
applications.  As  part  of  this  program,  TSC  is  conducting  ana- 
lytical and  experimental  studies  directed  towards  noise  and 
vibration  reduction  in  urban  rail  systems. 
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The  UMTA  Rail  Noise  Abatement  effort  will  bring  to- 
gether and  improve  existing  and  new  elements  into  a unified 
technology  consisting  of:  design  criteria  for  establishing 

goals,  noise  control  theory,  design  methods,  test  procedures, 
and  appropriate  documentation.  The  program  has  been  organized 
into  four  concurrent  and  interrelated  parts  which  will  be 
closely  coordinated  with  each  other  by  TSC.  They  are; 

1.  Assessment  of  Urban  Rail  Noise  and  Vibration 
Climates  and  Abatement  Options; 

2.  Test  and  Evaluation  of  State-of-the-Art 
Urban  Rail  Noise  Control  Techniques; 

3.  Wheel/Rail  Noise  and  Vibration  Control 
Technology ; 

4.  Track  and  Elevated  Structure  Noise  and 
Vibration  Control  Technology. 

This  report  is  a result  of  work  performed  under  part  4 of  the 
program  and  deals  with  the  development  and  engineering  use  of 
a model  for  predicting  the  effectiveness  of  floating-slab  track 
in  reducing  the  vibration  and  noise  transmitted  to  buildings 
near  subway  tunnels.  Results  of  field  measurements  for  a length 
of  floating-slab  track  on  the  New  York  City  transit  system  are 
also  presented. 

1.2  Program  Objectives 

In  an  effort  to  reduce  the  vibration  transmitted  from 
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subway  tunnels,  a number  of  techniques  have  been  suggested  and 
tried  [2] . The  most  successful  of  these  has  been  to  isolate 
the  vibrations  of  the  rail  from  the  tunnel  floor  by  means  of 
resilient  rail  fasteners  and  floating-slab  track.  Although 
the  floating-slab  tracks  appear  to  be  more  effective  than 
resilient  rail  fasteners  in  preventing  the  transmission  of 
vibration,  past  installations  have  been  so  expensive  that 
their  use  may  not  be  cost  effective. 

The  goal  of  the  work  presented  in  this  report  is  to 
develop  a model  for  predicting  the  performance  of  floating- 
slab  track.  It  is  hoped  that  by  developing  a basic  understanding 
of  these  devices,  the  transit  industry  can  implement  less  ex- 
pensive and  equally  effective  designs. 

1,3  State-of- the-Art  in  Floating-Slab  Track  Design 

Within  the  present  state-of-the-art  floating-slab 
tracks  are  designed  according  to  a simple  single-degree-of- 
freedom  model  shown  in  Fig.  1.1.  The  resonance  frequency  of 
the  model  is  set  to  equal  the  rigid  body  resonance  of  the 
floating  slab  on  its  supports. 

The  damping  is  determined  from  the  damping  properties 
of  the  slab  supports.  As  seen  in  Fig.  1.2  this  simple  model 
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FORCE  ON  RAIL  ( q) 


SLAB  MASS 
PER  UNIT 
AREA(yO^) 

T" 


DAMPING  ta 


I 


SLAB  SUPPORT 
STIFFNESS  PER 
UNIT  AREA  (Kg) 


FORCE  ON  THE 
TUNNEL  FLOOR 


b)  SINGLE  DEGREE  OF  FREEDOM  SYSTEM 


FIG.  1.1  SIMPLE  MODEL  USED  TO  PREDICT 
FLOATING  SLAB  VIBRATION  REDUCTION 
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PARAMETERS 

Tp  FORCE  TRANSMISSIBLITY,  RATIO  OF  THE  AMPLITUDES 
OF  THE  FORCE  ON  THE  RAIL  AND  THE  FORCE 
ON  THE  TUNNEL  FLOOR; 

f^s  RESONANCE  FREQUENCY  OF  THE  SLAB  ON  ITS 
SUPPORTS,  fos^=  K,/p,\ 

71  VISCOUS  DAMPING  LOSS  FACTOR  ( FOR  SOLID  TYPE 
''  DAMPING,  Ss  , LET  fog/f); 


7j^  = 2^  WHERE  ^ IS  100  TIMES  THE  PERCENT  OF 
CRITICAL  DAMPING  VALUE 


FIG  1.2  FORCE  TRANSMISSIBILITY  FOR  THE 
SINGLE  DEGREE  OF  FREEDOM  SYSTEM 
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predicts  that  the  forces  applied  to  the  tunnel  floor  are  ef- 
fectively reduced  at  frequencies  well  above  the  resonance 
frequency  and  are  increased  at  frequencies  near  the  resonance 
frequency.  A design  tradeoff  must  be  made  in  selecting  the 
damping.  High  damping  reduces  the  force  amplification  at  the 
resonance  frequency  but  at  the  expense  of  higher  frequency 
performance . 

Most  floating  slab  designs  have  set  the  resonant 
frequency  in  the  range  10  to  30  Hz.  The  objective  has  been  to 
set  the  resonance  frequency  as  low  as  possible  within  material 
and  static  displacement  constraints.  One  problem  has  been  to 
find  resilient  materials  that  can  withstand  the  high  static 
loads  and  severe  environmental  conditions  associated  with  tran- 
sit systems.  Most  designs  make  it  difficult  to  replace  the 
slab  supports  so  that  their  life  must  be  very  long. 

The  single-degree-of-f reedom  model  does  not  take  into 
account  the  flexibility  of  the  slab  or  effects  of  using  a resil- 
ient rail  fastener  between  the  rail  and  the  slab.  In  spite  of 
this  neglect  the  performance  predicted  by  the  simple  model  for 
a floating  slab  track  in  New  York  City  agrees  reasonably  well 
with  measured  data  that  we  have  taken.  We  believe  this  to  be  a 
fortuitous  result  for  the  particuJar  slab  design  that  was  studied. 
A more  realistic  model  may  be  required  for  other  designs. 
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In  this  report  we  develop  a new  model  that  takes 
into  account  the  flexibility  of  the  slab  and  the  effects  of 
resilient  rail  fasteners.  The  model  can  be  used  to  study 
both  long  continuous  slabs  and  resiliently  supported  ties. 

1.4  Summary  of  Results 

The  prediction  model  developed  in  this  report  allows 
the  determination  of  the  reduction  in  the  dynamic  forces  applied 
to  the  tunnel  floor.  A method  to  calculate  the  response  of  the 
tunnel  to  the  applied  forces  has  not  been  found  and  should  be 
the  subject  of  future  work.  Similarly,  methods  to  calculate 
the  transmission  of  tunnel  vibration  through  the  ground  to 
neighboring  buildings  have  not  been  developed  as  part  of  the 
work  carried  out  in  this  report.  However,  in  spite  of  these 
limitations  the  authors  believe  that  the  reduction  in  the  dynamic 
forces  acting  on  the  tunnel  floor  is  a good  measure  of  performance 
and  can  be  used  to  compare  different  floating-slab  designs. 

The  prediction  model  indicates  that  the  ratio  of  slab 
support  stiffness  to  slab  mass  should  be  as  small  as  possible 
within  imposed  limitations  on  maximum  allowable  deflections  during 
the  passby  of  a fully  loaded  train.  The  slab  performance  at  the 
rigid  body  resonance  of  the  slab  on  its  supports  is  governed  by 
the  damping  of  the  supports,  with  lightly-damped  supports  giving 
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poor  performance.  However,  a tradeoff  exists  in  selecting  the 
support  damping,  since  high  values  of  damping  cause  the  performance 
above  the  rigid  body  resonance  to  decrease. 

The  role  of  slab  support  damping  at  frequencies  above 
the  rigid  body  slab  resonance  is  varied.  For  very  low  values  of 
damping  performance  of  the  slab  is  poor  because  of  the  excitation 
of  lightly  damped  bending  waves  in  the  slab.  For  somewhat  higher 
values  of  support  damping  the  bending  wave  motion  of  the  slab  is 
reduced  to  a minimum  value  determined  largely  by  the  slab  mass. 
Finally,  for  large  values  of  damping  the  forces  transmitted  to 
the  tunnel  floor  increase  as  shown  by  the  simple  prediction 
chart  in  Fig.  1.2 

At  high  frequencies,  above  the  rigid  body  resonance 
of  the  rail  on  its  fasteners,  it  is  possible  for  the  free  bending 
wavespeed  in  the  rail  and  the  slab  to  coincide.  This  coincidence 
effect  reduces  the  performance  of  the  slab  at  these  frequencies 
and  should  be  avoided. 

The  field  study  results  support  the  general  conclusions 
of  the  prediction  model.  They  show  that  the  floating  slab  has 
a detrimental  effect  at  and  near  the  rigid  body  slab  resonance. 

For  the  particular  system  studied  the  vibration  levels  of  the 
tunnel  are  increased  by  the  floating  slab  approximately  10  dB 
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in  the  one-third  octave  band  centered  on  the  rigid  body  reso- 
nance frequency.  Above  this  resonance  frequency  the  performance 
of  the  slab  improves  rapidly  with  increasing  frequency.  For 
this  particular  slab,  the  support  damping  appears  to  be  suffi- 
ciently high  that  the  vibration  due  to  free  bending  waves  in 
the  slab  is  not  important  in  determining  the  slab  performance. 

We  conclude  from  this  study  that  floating  ties  can  be 
just  as  effective  as  long  floating  slabs.  The  floating  ties 
prevent  the  occurrence  of  coincidence  effects  and  eliminate 
free  bending  wave  response  of  the  track  support — both  beneficial 
effects  with  regard  to  reducing  the  dynamic  forces  transmitted 
to  the  tunnel  floor.  On  the  other  hand,  the  floating  ties  allow 
the  static  deflection  during  a train  passby  to  be  larger  and 
thereby  increase  the  maximum  stresses  on  the  rail. 

1.5  Organization  of  the  Report 

The  report  is  organized  into  three  basic  sections.  In 
Section  2 we  describe  the  mathematical  formulation  of  the  pre- 
diction model.  The  basic  approach  is  presented  and  predictions 
applicable  to  a high  and  low  frequency  regime  are  developed.  In 
Section  3 we  present  the  results  of  a field  study  of  the  performance 
of  a length  of  floating-slab  track  on  the  New  York  City  transit 
property.  Finally,  in  Section  4,  we  summarize  the  results  and 
show  how  they  can  be  applied  to  the  practical  engineering  design 
of  floating-slab  track. 
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2.  MATHEMATICAL  FORMULATION 


2 . 1 Introducti on 

This  section  presents  the  mathematical  formulation 
of  a model  to  predict  the  dynamics  of  floating-slab  track  systems 
and  the  influence  of  design  parameters  on  their  effectiveness  in 
reducing  vibration  transmission  from  trains  to  subway  floors  and 
hence  to  neighborhood  buildings. 

The  rationale  for  the  "floated  track  slab"  is  similar 
in  philosophy  to  a vibration  isolation  method  used  in  industrial 
situations  which  involves  mounting  vibrating  machinery  on  an 
inertia  block,  a reinforced  concrete  block  resiliently  mounted 
on  a floor.  In  a recent  study  [3] , Bender  summarizes  the  signi- 
ficant differences  between  the  floated  slab  system  and  a conven- 
tional inertia  block.  However,  his  work  is  principally  directed 
toward  the  investigation  of  acoustic  radiation  from  the  rail-slab 
system.  This  is  accomplished  by  a detailed  consideration  of  the 
separate  behaviors  of  rail  and  slab  and  a qualitative  analysis 
of  the  effects  of  dynamic  coupling  by  assuming  the  various  system 
elements  to  be  cascaded. 

The  present  effort  is  primarily  devoted  to  an  under- 
standing of  the  vibration  transmission  properties  of  the  floated 
slab  system.  In  addition,  the  results  of  the  previous  study  [3] 
have  been  extended  by  removing  certain  approximations  and 
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restrictive  assumptions.  More  specifically,  our  effort  provides 
a detailed  examination  of  the  dynamic  coupling  between  rail  and 
slab  and  an  approximate  solution  to  the  fully  coupled  equations 
of  motion  including  a Bernoulli-Euler  plate  formulation  for  the 
slab  and  visco-elastic  damping  associated  with  rail  fasteners 
and  slab  supports.  The  results  are  also  shown  to  apply  to 
floated  tie  systems.  We  shall  evaluate  appropriate  measures 
of  the  force  transmitted  to  the  floor  of  the  subway  tunnel  and 
the  acoustic  radiation  from  the  slab.  Guidelines  for  optimizing 
the  vibration  isolation  performance  are  presented. 

2.2  Basic  Formulation 

The  floated  slab  system  consists  of  rails  resiliently 
fastened  to  a concrete  slab  which  is  in  turn  mounted  on  the  tun- 
nel floor  by  resilient  pads.  The  basic  physical  model  to  be 
adopted  for  this  analysis  along  with  definitions  of  the  perti- 
nent parameters  are  presented  in  Fig.  2.1.  Representative  values 
of  these  parameters  are  given  in  Appendix  A. 


For  the  purpose  of  developing  guidelines  for  floated 
slab  design  it  is  deemed  sufficient  to  examine  the  interaction  of 
the  slab  with  a single  rail  excited  by  the  motion  of  one  wheel  of 
the  rail  vehicle.  To  provide  a basis  for  more  refined  analysis 
it  is  assumed  initially  that  both  rail  and  slab  are  infinite 


n 


= rail  mass  per  unit  length 

E = modulus  of  elasticity  of  the  rail 
material 

= moment  of  inertia  of  the  rail 

= slab  mass  per  unit  area 

E = modulus  of  elasticity  for  the  slab 
material 

I = moment  of  inertia  for  the  slab 
s 

D = slab  flexural  rigidity  per  unit  width 

h = slab  thickness 

L = slab  width 

y^  = position  of  rail  fasteners  on  the  slab 

q = applied  force  due  to  rail-rail  inter- 
action 

* 

K = visco-elastic  stiffness  modulus  for  rail 
fasteners  per  unit  length  of  rail 

* 

K = visco-elastic  stiffness  modulus  of  slab 

3 

supports  per  unit  area  of  slab 
For  harmonic  response; 

K*=K  (l-j6),  K*=K  (l-j6) 
r r j > g g ' s 

where  6^  and  6^  are  the  solid  type  damping  loss 
factors  and  j = /-I . 


FIG.  2.1  SLAB  GEOMETRY  AND  DEFINITIONS  OF  PARAMETERS 


12 


along  the  track  and  that  the  slab  is  of  constant  rectangular 
cross-section  and  homogeneous  material  properties.  In  Section 
2.10  the  analysis  is  extended  so  that  it  can  be  applied  to  floating- 
tie  tracks.  Also,  only  the  vertical  motions  of  rail  and  slab  are 
considered.  Damping  is  accounted  for  in  the  formulation  by  ideal- 
izing the  slab  mounts  and  rail  fasteners  by  continuous  visco- 
elastic supports. 

In  all  practical  instances,  the  train  speed  is  much  smaller 
than  the  wave  speed  of  vibrations  transmitted  down  the  rail  and 
slab.  Thus,  the  force  due  to  wheel /rail  interaction  may  be 
modeled  by  a stationary,  vertical  point  force  located  at  the  origin 
of  the  X coordinate  [4] . 

The  fluid  loading  of  rail  and  slab  due  to  the  acoustic  space 
is  neglected  so  that  system  dynamics  are  calculated  in  vacuo.  In 
some  cases  of  practical  interest,  however,  the  effective  stiffness 
of  air  enclosed  under  the  slab  must  be  included  in  calculating  the 
stiffness  of  the  slab  supports. 

Finally,  the  dynamic  interaction  of  the  slab  with  the  tunnel 
floor  is  neglected  in  this  analysis.  That  is,  the  tunnel  floor  is 
taken  to  be  a rigid  support.  This  assumption  is  rather  restrictive 
as  it  requires  the  tunnel  floor  impedance  to  be  much  larger  over 
all  frequencies  than  the  impedances  of  all  other  system  components. 
However,  it  is  anticipated  that  once  a determination  of  actual  tun- 
nel floor  impedance  has  been  made,  the  present  results  can  be  suit- 
ably extended  so  that  the  tunnel  floor  vibration  can  be  predicted. 
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The  specific  analytical  idealizations  used  here  in  order  to 


make  the  analysis  tractable  are 
System  Element 

Rail 

Slab 

Rail  fasteners  and  slab  mounts.. 


Wheel-Rail  interaction 


summarized  below. 

Analytical  Model 

Bernoulli-Euler  beam 
Bernoulli-Euler  plate  with 
edges  Y=0,L  assumed  stress 
free . 

Continuous,  viscoelastic 
supports  with  frequency 
independent  loss  factors,  i.e., 
for  harmonic  response  these 
are  modeled  with  complex 
elastic  moduli. 

Point  force  on  the  rail  at 
x=0.  The  force  time  histories 
are  assumed  stationary  random 
processes . 


Results  are  derived  for  two  kinds  of  power  spectral  den- 
sities of  the  existing  excitation  at  the  wheel/rail  contact  points 
(i)  a flat  spectral  density  and  (ii)  the  spectrum  arising  from  the 
micro-roughness  model  of  wheel-rail  interaction  [5] . 
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Before  presenting  the  picture  of  rail  and  slab  response 
developed  under  the  above  assumptions,  we  shall  discuss  the 
dynamic  characteristics  of  various  individual  system  elements 
and  outline  our  general  approach. 

2.3  Characteristics  of  Individual  Rail  and  Slab  Motions 

Some  appreciation  of  the  overall  response  characteristics 
can  be  obtained  by  considering  the  unforced  and  undamped  motion  of 

(1)  the  rail  on  its  fasteners  with  the  slab  treated  as  rigid  and 

(2)  the  slab  on  its  mounts  without  the  rail  and  fasteners.  We 
shall  call  those  two  systems  the  uncoupled  rail  and  uncoupled  slab , 
respectively . 


First,  consider  the  free,  undamped  motion  of  the  un- 
coupled rail.  The  rail  displacement  in  this  case,  5(t,x),  is  in 
the  form  of  sinusoidal  bending  waves. 


C(t,x) 


j (k  x±o)t) 
e X 


(2.1) 


where  the  wave  number,  k^,  is  related  to  radian  frequency,  oj, 
through  the  dispersion  relation. 


(2.2a) 


15 


where  the  wavenumber  for  freely  propagating  waves  in  the  rail, 

k , is 
r 


r 


[()^)^-l] 

(JL) 

or 


(2.2b) 


and  (jj  , the  natural  frequency  of  the  rail  moving  as  a rigid 
or 

body  on  its  fasteners,  is  given  by 


ixi 


or 


/K  /p 
r 


r 


and  the  parameter  is  given  by 


(2 . 2c) 


1 El 
r r 


(2 . 2d) 


where  is  the  effective  rail  fastener  stiffness  per  unit  length 

of  rail  and  El  is  the  bending  stiffness  of  the  rail.  This  dis- 
r r 

persion  relation  is  shown  in  Fig.  2.2  for  three  different  values 
of  fastener  stiffness.  Clearly  the  character  of  the  rail  motion 
depends  on  the  value  of  oo  relative  to  w . For  m > w , free  bend- 
ing  waves  can  propagate  down  the  rail,  while  for  w < propa- 
gating bending  waves  cannot  exist.  Finally  for  co  = the  rail 

oscillates  as  a rigid  body  on  its  fasteners. 


The  situation  in  the  case  of  the  slab  is  more  complicated. 
The  free  motion  of  the  slab  w(t,x,y)  may  be  expressed  as  a combina- 
tion of  motions  of  the  form 
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FREQUENCY  IN  HERTZ 


0 2 4 6 8 10 


DIMENSIONLESS  WAVENUMBER  kr  L 


FIG.  2.2  RAIL  DISPERSION  CURVES  FOR 
THREE  VALUES  OF  FASTENER  STIFFNESS 
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(2.3) 


where  the  (n=l , 2 , . . . ) are  the  cross-wise  normal  modes  of  the 

slab  displacement.  Each  mode  is  seen  to  propagate  along  the 
slab  as  a sinusoidal  wave.  These  mode  shapes  depend  upon  the 
wavenumber  k and  the  Poisson  ratio,  u,  as  well  as  upon  y. 
However,  for  y small,  as  is  typical  of  concrete  (y  ==  .15)  , the 
ijj^'s  are  approximately  independent  of  y. 


number  and  frequency  must  be  related  by  a distinct  dispersion 
relation.  For  the  values  of  parameters  given  in  Appendix  A, 
they  are  given  to  within  2%  accuracy  by: 


Furthermore,  for  each  slab  normal  mode,  the  wave- 


Y 


n 


(2.4) 


where  a is  related  to  the  frequency  by 


a 


4 


1) 


(2.5a) 


os 
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where  the  natural  frequency  of  the  slab  moving  as  a rigid 

body  on  its  resilient  pads,  is  given  by 


(2.5b) 


and  the  parameter  is  given  by 


K 


2 


or 

CO 


os 


(2.5c) 


where  is  the  effective  stiffness  of  the  slab  supports  per  unit 
area  of  the  slab  and  where  D is  the  bending  stiffness  of  the  slab, 
per  unit  width,  D=  E^h^/12,* 


Equations  2.4  are  shown  in  Fig.  2.3  (using  the  values 
given  in  Appendix  A)  along  with  the  rail  dispersion  curves  given 
earlier.  As  in  the  case  of  the  rail,  there  is  a cutoff  frequency 
for  each  slab  mode,  below  which  that  mode  cannot  propagate  freely 
down  the  slab.  These  cutoff  frequencies,  (n=l,2...)  are 

given  approximately  by 


ns 


os 


= 1 + 


(0 

or 

2 

(0 

os 


n = 1 


1 1 

— (y(n-  i)) 
<2  L 2 


n > 1 


(2.6) 


*In  this  relation  we  have  neglected  the  small  effect  of  Poisson 
coupling  which  acts  to  increase  the  bending  stiffness  in  a plate 
2 -1 

by  a factor  (1-y  ) where  y is  Poisson's  ratio. 
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FIG.  2.3  DISPERSION  CURVES  FOR  THE  FIRST 
FOUR  CROSS  MODES  OF  A CONCRETE  SLAB 
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gives  the  cross-sectional  shape  of  the  dis- 
~th 

placement  associated  with  the  n mode  and  is  a function  of  k^. 

The  first  mode  shape,  merely  represents  a rigid  body 

translation  of  the  slab  cross-section.  As  shown  in  Figs.  2.4 

and  2.5  the  higher  mode  shapes  depend  upon  k^.  These  plots 

show  the  second  and  third  mode  shapes  for  very  low  to  very 

high  wavenumbers.  These  functions  have  been  normalized  to  -1 

at  the  slab  free  edge.  It  is  seen  that  for  very  low  wave- 

th 

numbers,  ij;^(n>l)  is  approximately  the  n mode  shape  of  a free- 
free  beam  on  an  elastic  foundation,  while  for  very  large  wave- 
numbers,  approaches  the  limit 

iL  - - cos  nir  n = 2,  3,...  (2.7) 

n L 

Contrary  to  what  has  been  stated  in  Ref.  [3]  the  slab 
modes  are  orthogonal.  That  is: 


,L 

0 


(y,k  ) \p  (y,k  )dy 
n X ^m  ^ X 


0 


(2.8a) 


for  n m.  This  is  demonstrated  in  Appendix  B.  Furthermore 
the  modes  can  be  so  defined  that  their  norm  is  L.  Thus,  we  may 
write 
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MODE  SHAPE  AMPLITUDE  ^ 


FIG.  2.4  SECOND  SLAB  MODE  SHAPE 
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MODE  SHAPE  AMPLITUDE 


I 


FIG.  2.5  THIRD  SLAB  MODE  SHAPE 
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rL 


(2.8b) 


il;  ij;  dy 


L , n=m 

'0, 


The  slab  mode  shapes  defined  above  may  be  used  to 
synthesize  the  actual  forced  and  damped  motion  of  the  slab, 
since  in  general  the  slab  displacement  may  be  expressed  as  a 
combination  of  functions  multiplied  by  appropriate  weighting 
factors . 


The  dispersion  curves  of  uncoupled  rail  and  slab 
can  give  a qualitative  picture  of  the  actual  response.  In  par- 
ticular, for  cases  wherein  the  dispersion  curves  of  rail  and 
slab  neither  cross  nor  approach  one  another  very  closely,  the 
system  elements  can  be  cascaded.  That  is,  the  forced  rail 
motion  is  calculated  ass\iming  the  slab  rigid  and  the  slab  motion 
is  taken  to  be  the  response  of  the  uncoupled  slab  to  forcing  by 
the  rail  motion  through  the  fastener  stiffness. 

When  rail  and  slab  dispersion  curves  cross,  however, 
a cascade  analysis  is  not  valid.  This  condition  is  called 
structural  wave  coincidence,  since  at  the  particular  frequency 
at  which  it  occurs,  the  wavelength  of  free  bending  waves  of 
both  rail  and  slab  are  the  same.  At  the  coincidence  frequency 
the  system  components  are  very  tightly  coupled  and  a detailed 
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analysis  of  the  coupling  dynamics  must  be  carried  out.  In  the 
following  we  outline  an  analysis  that  does  just  this. 


2.4  Rail-Slab  Interaction  and  Total  System  Response 

Now  consider  the  dynamics  of  the  slab  fully  coupled  to 
the  rail.  Our  basic  approach  is  as  follows. 

The  slab  displacement  can  be  expressed  as  a combination 
of  its  normal  modes: 


where  the  w are  the  modal  coordinates  of  the  slab  motion.  The 
n 

next  step  is  to  express  both  slab  and  rail  displacements  as  a 
combination  of  components,  each  of  which  represents  a displace- 
ment varying  sinusoidally  with  x at  a given  wavenumber.  The 
mathematical  statement  is: 


CO 


w(x,y,t)  = Z w^(x,t)  ’Jjjj(y) 


n=l 


(2.9) 


dk 


X 


(2.10) 


— 00 


dk  e 


jk  X 
X 


w(k^,y , t) 


(2.11) 


X 
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Or,  on  substituting  Eq.  2.9  into  Eq.  2.11,  we  get 


w(x,y,t)  = Z 


n=l 


2ir 


dk 


jk  X 
X 


X 


»„(kx,t)l 


(2.12) 


r\j 

g and  w are  seen  to  be  the  Fourier  transforms  of  C w wxth 

respect  to  x.  These  quantities  might  also  be  termed  the  wave- 
number  components  of  rail  and  slab  motion. 

Eq.  2.12  can  be  interpreted  as  first  a decomposition 
of  the  slab  displacement  into  its  normal  modes  and  then  a re- 
solution of  each  modal  coordinate  into  its  wavenumber  components. 

The  solution  procedure  is  to  obtain  ^(k^,t)  and  the 

w (k  ,t)  in  terms  of  the  exciting  force,  q(t) , and  then  to 
n ^ 

determine  ?(k  ,t)  and  w (x,t)  for  all  values  of  n by  calculating 
X n 

inverse  Fourier  transforms. 

Substitution  of  Eqs.  2.10  and  2.12  into  the  partial 
differential  equations  of  motion  results  in  an  infinite  set 
of  coupled  ordinary  differential  equations  for  ^(k^,t)  and 
w^(k^,t),  the  y dependence  being  completely  eliminated.  The 
substitution  procedure  is  discussed  in  Appendix  C.  It  comes 
about  that  for  a given  k^,  the  wavenumber  components 
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of  the  rail  and  slab  modes  behave  as  the  displacements  of  a 
system  of  coupled  mechanical  oscillators,  one  oscillator  being 
associated  with  each  slab  mode  and  with  the  rail  response. 

To  a good  approximation  these  oscillators  are  coupled  together 
only  by  spring  forces  arising  from  the  stiffness  of  rail  fasteners 
and  slab  supports. 

The  analogy  which  comes  out  of  our  analysis  is  shown 
in  Fig.  2.6  where  only  the  first  two  slab  modes  are  considered. 

For  all  k^,  the  rail  oscillator  is  excited  by  the  driving  force 
on  the  rail,  q(t),  while  the  slab  mode  oscillators  are  only  in- 
directly excited  by  their  coupling  to  the  rail.  As  the  figure 
illustrates,  each  slab  mode  oscillator  is  not  only  connected  to 
the  rail  oscillator,  but  is  coupled  to  all  the  other  slab 
oscillators . 

For  a given  k^,  the  resonant  frequencies  of  the  un- 
coupled oscillators  are  given  simply  by  the  dispersion  curves 
of  the  uncoupled  rail  and  slab.  It  should  be  noted  that  all 
parameters  shown  in  Fig.  2.6  are  functions  of  the  wavenumber 
and  are  always  different  for  different  wavenumbers.  Not  shown, 
but  implicit  in  the  figure,  is  the  frequency  independent  visco- 
elastic damping  associated  with  the  springs  K , K , K , etc. 

r s^  S2 
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RAIL 

OSCILLATOR 


V SLAB 

f oscillator; 


j 


All  masses  and  stiffnesses  are  functions  of  k : 

X 


= -T-'  ^'sl  = ^s2 

CO 

or 


k = k , k 
rr  X rs 


^sl 


^1  ,,  4 , “os  *, 


GO 


or 


4 4 

(k>Yo  + 


‘s2  . , ,1  , 2 3 '"'’x ' ' 2 


os  * , 
2 '^2^ 
^or 


where  for  harmonic  response 

K*  = - j«^)  , K*  = K^(l  - j6^). 

The  rigid  connecting  bar  is  free  to  rotate  about 
pivot  a and  to  translate  vertically. 


FIG.  2.6  fiECHANICAL  OSCILLATOR  ANALOGY  FOR  THE  FLOATED  SLAB  SYSTEM 
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First,  let  us  consider  the  free  motion  of  this  com- 
posite system.  It  is  found  that  the  free  motion  of  each  os- 


cillator is  a combination  of  the  normal  modes  of  the  total 
coupled  system  (i.e.,  those  motions  in  which  all  displacements 
vary  sinusoidally  with  the  same  frequency) . The  frequency  of  each 
mode,  when  considered  as  a function  of  yields  a distinct  dis- 
persion curve.  A typical  situation  is  shown  schematically  in 
Fig-  2.7,  which  shows  two  such  curves  along  with  the  uncoupled 
rail  and  slab  dispersion  curves.  For  k <<  k being  the 

coincidence  wavenumber)/  the  lower  frequency  mode  nearly  coincides 
with  the  uncoupled  slab  while  the  higher  mode  follows  the 

rail.  For  k >>  k the  correspondence  is  reversed.  At  w=a) 

X xc  c 

the  two  dispersion  curves  have  a minimum  separation  along  the 
k axis  equal  to  A , where 


where  and  are  defined  in  Eqs.  2. 2d  and  2.5c,  S is  the 
ratio  of  rail  mass  per  unit  length  to  the  slab  mass  per  unit 
length , 


X 


n 


(2.13a) 


s 


(2.13b) 


th. 

and  (y  ) is  the  value  of  the  n^  slab  mode  shape  at  the  rail 
n o 


29 


FREQUENCY 


FIG.  2.7  COUPLED  MODE  DISPERSION  CURVES 
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position.  The  minimum  separation , , is,  in  fact,  a measure  of 

th 

the  effects  of  coupling  forces  between  the  rail  and  the  n slab 
mode.  It  should  be  noted  that  since  we  have  lumped  all  damping 
into  the  coupling  stiffnesses  which  are  small  in  comparison  with 
bending  forces,  the  damping  band-width  of  each  oscillator  is 
always  smaller  than  the  coupling.  This  is  a basic  condition  for 
the  validity  of  many  of  the  results  presented  here. 

Returning  to  our  coupled  oscillator  analogy,  note  that 
at  wavenumber  coincidence  the  rail  and  slab  behave  as  identical, 
lightly  coupled  oscillators.  Thus,  since,  as  judged  from  the 
typical  values  given  in  Appendix  A , is  quite  small,  the 
overall  amplitude  of  forced  response  is  likely  to  be  relatively 
large  at  the  coincidence  frequency. 

Now  consider  the  response  to  a harmonic  point  force 
excitation  of  the  rail 

q(t)  = q (2.14) 

where  we  shall  take  q to  be  constant  over  frequency  during  this 
discussion.  The  resulting  wave  number  components  of  rail  and 
slab  modes  are  also  harmonic  with  the  same  frequency! 


(2.15a) 


w (k  ,t)  = w (k  ,o))e  n = 1,2,...  (2.1S.j) 

n X n X 


When  these  relations  are  placed  in  the  infinite  set 

of  coupled  oscillator  equations,  the  exact  solutions  for 

5 (k  ,co)  and  w (k  ,o))  can  be  determined.  Finally,  the  inverse 
X n X 

/s 

spatial  Fourier  transforms  of  C and  can  be  obtained  approxi- 
mately to  give  the  spatial  distribution  of  response 

C(x,t)  = C (x,(d)  (2.16a) 


'X, 

w (x  , t ) = w (x  , (jj)  e 
n n 


- j03t 


n = 1 , 2 , . . . 


(2.16b) 


% 

'X>  '\j  ^ ^ 

where  ? and  w are  the  spatial  Fourier  transforms  of  E and  w . 

n n 


The  spatial  distributions  of  rail  and  slab  motion  have 
been  obtained  in  this  way  under  the  ass\omption  that  only  the 
first  slab  mode  is  significantly  excited.  This  restriction 
leaves  out  the  effects  due  to  participation  of  higher  modes, 
but  nevertheless  adequately  illustrates  the  salient  aspects 
of  response. 
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are 


% 

% a. 

The  resulting  expressions  for  ^(x)  and  Wj^(x) 
extremely  lengthy  and  not  relevant  to  the  actual  design  of  floating 
track  slabs,  so  we  shall  merely  summarize  the  qualitative  results 
here.  The  expressions  are  given  in  Appendix  C. 

(a)  For  0)  bending  waves  propagate,  the 

displacements  of  rail  and  slab  consisting  of  near- 
field components  only.  For  to  <to  < oj  the  near 
field  of  the  rail  excites  the  far  field  of  the 

slab.  Above  to  and  to  , the  far  fields  of 
or  os 

both  rail  and  first  slab  mode  consist  of  two 
component  waves  propagating  at  the  rail  free 
bending  wave  speed  and  at  the  slab  wave  speed. 

As  long  as  to  is  not  within  a damping  bandwidth 
of  to  or  to  , the  role  of  damping  is  limited 
to  the  attenuation  of  bending  waves. 

(b)  For  to  = to^g  the  overall  amplitudes  of  both  rail 

and  slab  reach  a resonant  maximum,  limited  primarily 
by  the  damping  in  the  slab  supports  and  secondarily 
by  damping  in  the  rail  fasteners  and  the  slab  it- 
self. The  shape  of  the  response  consists  essentially 
of  a nearfield  with  maximum  value  at  the  driving 
point  extending  a large  distance  from  the  driving 
point.  Smaller  bending  wave  components,  propagating 
at  the  rail  wave  speed  also  exist.  For  o)  = oj  , 
these  remarks  also  apply  if  the  roles  of  rail  and  slab 
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are  interchanged. 


(c)  For  a structural  wave  coincidence  occurring  at 

the  driving  frequency  the  maximum  response  of 

the  rail  occurs  at  the  driving  point  while  the 

slab  response  is  very  small  there  (of  order 
1 I 4 

At/  k , where  k is  allowed  to  be  complex  to 
allow  for  damping  in  the  rail  fasteners.  For 
moderate  distances  from  the  driving  point,  the 
far  field  components  of  response  are  approximately: 


1 

4E  I 
r 


j k X 

e ' r ' 


cos 


j 


(2.17a) 


sin 


(2.17b) 


The  bending  waves  are  thus  modulated  sinusoidally, 

the  modulation  of  rail  and  slab  being  out  of  phase. 

This  is  essentially  a "beating"  phenomenon  where 

the  distance  between  successive  "beats"  is  ^lk  1^. 

A ' r ' 

Energy  is  continually  being  passed  from  rail  to  slab 
and  back  again  as  the  bending  waves  propagate  down 
the  track.  For  distances  sufficiently  large  that 
vibration  attenuation  due  to  damping  becomes  signifi- 
cant, Eqs.  2.17  are  no  longer  valid.  At  these 
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distances  "beating"  subsides,  the  vibration 


amplitudes  reach  a steady  state  and  rail  and 
slab  motions  eventually  come  to  be  in  phase. 

The  amplitudes  of  response  to  a harmonic  driving  force 
can  easily  be  used  to  obtain  the  driving  point  impedance  of  the 
rail,  z^(u)).  Again,  this  has  been  obtained  by  assuming  inter- 
action of  the  rail  with  the  first  slab  mode  only.  Neglecting 
terms  of  the  order  of  the  damping  gives: 


4|a|3  |k^|3 


(2.18) 


a|a|^+  b|k^l^ 


where 


k 


4 


< 0 


r 


4 4 2 

a = (k  - a + 

r 2 1 


(2.19a) 


/2,  < 0 


b = - X< 


> 


> 0 


J 


(2.19b) 


(2.19c) 


where  \p^  denotes  (y  ) . 

ri  ^1  -'o 
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As  can  be  seen  from  Fig.  2.8,  which  shows  |z^|  plotted 

versus  frequency  for  a typical  case,  the  magnitude  of  the  rail 

impedance  shows  a resonant  minimum  at  o)  as  well  as  at  w 
^ os  or 

On  the  whole,  is  very  close  to  the  driving  point  impedance 
of  the  uncoupled  rail,  except  within  a narrow  region  centered 
about  0)^^  where  its  magnitude  drops  off  sharply.  When  inter- 
action with  higher  slab  modes  is  accounted  for,  z^  is  found  to 
have  resonant  minima  at  each  slab  mode  cutoff  frequency  and  to 
be  close  to  the  uncoupled  rail  impedance  elsewhere. 

To  conclude  this  section  on  overall  system  response 
and  the  solution  procedure,  consider  the  response  to  a stationary 
random  excitation.  In  this  case  response  spectra  may  be  immedi- 
ately obtained  in  terms  of  the  harmonic  response  and  the  excita- 
tion spectrum.  For  example,  the  power  spectral  densities  of  rail 
and  slab  displacements  at  fixed  locations  are: 


^ 2 

S (w;x)  = S„(w) 

q ^ 


“ w (x,w)  2 

S^(w;x,y)  = \l  ip^(y)  -IL_ 1 


n=l 


(2.20) 


(2.21) 


where  and  are  the  spectral  densities  of  rail  response  at 
X and  the  slab  response  at  x and  y,  respectively,  and  is  the 
spectral  density  of  the  driving  force. 
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NORMALIZED  IMPEDANCE 


FIG.  2.8  DRIVING  POINT  IMPEDANCE  OF  THIS  RAIL 
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Although  the  spectra  of  the  rail  and  slab  responses 
at  a particular  point  are  of  some  interest  academically , they 
do  not  form  a convenient  basis  for  engineering  design  and  will 
not  be  calculated  in  this  report.  Instead  we  will  develop  and 
use  a measure  of  the  total  dynamic  force  acting  on  the  tunnel  floor. 


2.5  Force  Transmitted  to  the  Tunnel  Floor 

Of  major  concern  is  the  extent  to  which  the  floated 
slab  system  isolates  the  subway  tunnel  from  the  dynamic  forces 
imparted  to  the  rail  due  to  passage  of  a rail  vehicle.  The 
pressure  (force  per  unit  area)  on  the  tunnel  floor  at  location 
(x,y)  is  related  to  the  displacement  of  the  slab,  w(x,y,t)  and 
the  effective  spring  constant  of  the  slab  supports. 

F(x,y,t)  = KgW(x,y,t)  (2.22) 

Of  course,  the  ultimate  question  concerns  the  vibration 
of  the  tunnel  floor  produced  by  F(x,y,t).  Although  the  response 
characteristics  of  the  subway  tunnel  are  as  yet  unknown,  it  is 
still  clear  that,  to  some  approximation,  the  vibration  isolation 
effectiveness  of  floated  slab  systems  can  be  optimized  by  making 
a suitable  measure  of  transmitted  pressure  as  small  as  possible. 

The  effectiveness  of  conventional  vibration  isolation 
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is  usually  measured  by  the  force  transmissibility  which  is  the 


magnitude  of  the  ratio  of  force  transmitted  to  a floor  or  sup- 
port to  a sinusoidal  force  on  the  isolation  system.  A similar 
measure  of  vibration  isolation  effectiveness  might  be  defined 
for  the  floated  slab. 

One  possible  definition  of  force  transmissibility  is 
the  ratio  of  the  force  per  unit  area  on  the  tunnel  floor  under 
the  driving  point  to  the  force  on  the  rail.  On  closer  inspection, 
however,  this  proves  inadequate  since  it  fails  to  account  for  the 
effects  of  the  spatially  distributed  nature  of  the  slab  response. 

For  example,  there  can  be  different  frequencies  at  which  the  slab 
response  at  the  driving  point  is  the  same  but  at  one  frequency 
the  response  is  limited  to  a small  area  near  the  driving  point, 
whereas  at  the  other,  the  response  extends  over  a very  large 
region.  Also,  as  noted  in  the  last  section,  the  slab  response  at 
the  driving  point  is  nearly  zero  at  coincidence.  Thus,  this  mea- 
sure of  transmissibility  would  be  very  small  even  though  the  overall 
force  on  the  tunnel  is  significant  over  a great  distance.  The 
same  remarks  might  be  made  concerning  any  similar  measure  based 
on  the  tunnel  floor  force  at  only  one  location  on  the  slab. 

A better  but  still  inconvenient  alternative  can  be  ob- 
tained by  considering  the  mechanical  oscillator  analogy  presented 
above.  Just  as  was  done  for  the  slab  displacement,  the  forces  on 
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the  tunnel  floor  may  be  separated  into  the  portions  attributable 
to  each  slab  mode  and  each  of  these  may  be  decomposed  into  their 
wavenumber  components,  i.e.,  for  a sinusoidal  force  on  the  rail: 


F(x,y,t)  = Z 2^ 
n=l  _ 


dk 


jk  X 

. X 


*„(y)F^(k^,a.) 


-j  0)t 

i ’■* 


(2.23) 


til 

where  F is  the  contribution  from  the  k wavenumber  component 
th 

of  the  n^  slab  mode.  Placing  Eq . 2.22  into  Eq.  2.23  and  com- 
paring with  Eq.  2.12,  it  is  seen  that 


/s 


(k^,m) 


/s 

= (l-j6^) w (k  ,w) 

s s n X 


(2.24) 


where  the  slab  support  damping  factor,  <5^,  has  been  included  by 
allowing  the  stiffness  to  be  complex,  = K^(l-j6^).  Eq.  2.24 
together  with  Eq.  2.23  relates  the  tunnel  floor  force  directly  to 
the  coupled  oscillator  analogy.  Thus,  a transmissibility  can  be 
defined  for  each  slab  mode  oscillator  at  a given  wavenumber  by 

F 

^ (2.25) 

q 

and  the  effectiveness  of  the  slab  might  be  gauged  by  considering 
this  function  for  each  slab  mode  and  wavenumber. 


X 


Fig.  2.9  shows  the  transmissibility  for  the  first  slab 
mode  at  zero  wavenumber  for  rails  with  a range  of  fastener  stiff- 
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DIMENSIONLESS  FREQUENCY  f/  for 


FIG.  2.9  CONSTANT  WAVENUMBER  TRANSMISSIBILITY 
CURVES  FOR  THE  FIRST  SLAB  MODE 
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i! 

f ' 

, i 

i 

i 

nesses.  The  first  slab  mode  shape  is  just  a constant  over  the  | 

cross-section  and  the  zero  wavenumber  component  is  that  part  of  ^ 

the  slab  displacement  consisting  of  a rigid  body  motion  all  along 

its  length.  Hence  this  plot  shows  a transmissibility  arrived  at 

by  conventional  vibration  isolation  theory  modeling  the  slab  as  j 

a rigid  mass.  Moreover,  the  result  indicated  by  these  plots  -- 

that  force  transmission  (above  slab  resonant  frequency)  decreases 

with  increasing  fastener  stiffness  is  just  what  would  be  obtained  ! 

with  a conventional  isolation  system.  Thus,  if  the  coupled  oscil-  , 

)' 

lators  of  Fig.  2.6  behaved  in  the  same  fashion  for  all  wavenumbers,'? 

I 

the  situation  would  be  simple  indeed.  Unfortunately  this  is  not 
the  case,  as  can  be  seen  from  Fig.  2.10.  The  figure  shows  Eq.  | 

2.25  for  the  second  slab  mode  and  the  intermediate  rail  dispersion  j 
curve  of  Fig.  2.3  plotted  for  various  wavenumbers.  As  the  wave- 
number  increases  from  zero,  the  distance  between  peaks  is  completely 
altered.  Moreover,  for  k L=4  (at  v/hich  a coincidence  occurs)  the 

X I [ 

peak  at  - — =;  2 0 is  some  5 dB  above  what  it  is  for  k L=0.  Eq.  2.25  ' 
os 

thus  gives  different  measures  of  transmissibility  for  each  wavenumbe^ 
component  of  force.  Hence,  it  can  be  concluded  that  the  above  mea- 
sure of  slab  effectiveness  is  also  inconvenient.  What  is  needed  is 
an  overall  measure  of  the  tunnel  floor  force.  To  define  such  a 

quantity,  we  begin  by  considering  the  case  of  random  excitation. 

. 'i 

Clearly  in  this  case  the  spectrum  of  the  tunnel  floor  force,  " 

S„(s,y,03),  is  a valid  measure  of  the  force  at  each  location  on  the 

r 

slab.  But  S (x,y,(jo)  varies  considerably  with  x and  y and  so  does 

F I 
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not  provide  a single  measure  of  slab  performance  as  a whole  at 


a given  o). 


An  overall  measure  of  force  on  the  tunnel  is  the 
integral  of  over  the  whole  slab  area,  i.e., 


np(o))  = 


J -c 


dy 


Sp  (s,y,o)) 


(2.26) 


We  shall  call  the  integrated  spectrum  of  the  tunnel  floor 

r 

force.  Using  Eqs.  2.24  and  2.21  and  the  orthogonality  relations 
for  the  ip  , we  obtain: 


np(m)  = S (OJ) 

n=l  ^ n 


% 

'\y 


n„  (to)  = K^(l  + 6^)  L S^(to) 
s s q 


r " 

W (x,(jO) 

dx 

' — 00 

n 

/N 

rr 

(2.27) 


IIp  has  the  advantage  that  it  includes  the  effects  of  far  field 
slab  response,  i.e.,  it  is  large  when  the  region  over  which  the 
slab  response  is  significant  becomes  large. 


By  Parsevals ' theorem  the  integration  over  x can  be 
changed  to  one  over  k . Thus, 


"f  <“)  = ^2(1+62)  L S (a.)  i 

n ^ y 


r 


dk 


w (k  ,(o) 
n X 


(2.28) 


2tt 


dk  S (k  ,(o) 

X ' X 
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oscillator 


t h 

where  S_  (k  ,o))  is  the  spectral  density  of  the  n 
F X 
n 

response  at  wavenumber  k . TI  is  thus,  in  a sense,  the  average 

X t 

of  all  such  spectra  over  all  wavenumbers.  Note  that  when  n„  is 

r 

divided  by  it  gives  a measure  of  overall  dynamic  force 

transmissibility . 

Finally,  it  will  be  seen  that  the  integrated  spectrum  of 
the  force  on  the  tunnel  floor  is  closely  related  to  an  approxima- 
tion for  the  total  acoustic  power  radiated  by  the  slab. 


2.6  Calculation  of  the  Integrated  Spectrum 

The  procedure  for  calculating  the  integrated  spectrum 
of  the  force  per  unit  area  acting  on  the  tunnel  floor  depends  on 
the  values  of  the  damping  loss  factors  for  the  slab  and  its  sup- 
ports. For  small  values  of  damping  we  need  consider  only  the 
resonant  response  of  the  slab  modes  of  vibration.  Calculations 
for  this  case  are  presented  in  Section  2.6.1.  For  larger  values 
of  damping  contributions  from  non-resonant  slab  response  must  also 
be  included  in  calculation  of  the  integrated  spectrum.  The  cal- 
culations for  large  damping  are  presented  in  Section  2.6.2. 


2.6.1  Light  Damping 


The  integrated  spectrum  has  been  calculated  from  Eq.  2.28 
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under  the  assumption  that  the  damping  loss  factors  of  each  con- 


stant k oscillator  are  much  smaller  than  the  coupling  loss 

X 

factors  between  oscillators.  It  is  worthwhile  to  discuss  at  some 


length 

the  approximation  used  in  the  calculation.  The  integrand 

w (k  ,m)  ^ 

1 

q 

is  a positive  function  over  the  whole  m plane.  For  light 

damping  this  function  has  large  maxima  near  the  dispersion  curves 
of  the  coupled  system  normal  modes.  Thus,  the  integrand  can  be 
thought  of  as  a distribution  of  mass  in  the  k^  - oj  plane,  with 
large  maxima  near  the  dispersion  curves  of  the  coupled  system 


normal 

modes.  If  damping  is  small  then  most  of  the  contribution 

to  the 

integral  comes  from  the  regions  near  the  modal  dispersion 

curves , 

. The  subsequent  procedure  for  calculating  the  lightly 

daraped 

response  is  illustrated  in  Fig.  2.11.  First  the  "spectral 

mass"  in  the  interval  [k  , k + dk  ] attributable  to  the  peaks  at 
each  dispersion  curve  is  calculated.  This  can  be  done  by  rela- 
tively simple  integrations  over  w.  Next,  the  spectral  masses  are 
lumped  onto  the  corresponding  dispersion  curves,  i.e.,  the  inte- 
grand is  treated  as  a set  of  line  densities  concentrated  solely  on 
the  dispersion  curves.  Finally,  these  spectral  masses  are  placed 
within  the  intervals  dw^  and  dm2  in  the  manner  indicated  in  the 
figure.  The  result  is  a spectral  distribution  which  is  correct 
only  within  frequency  intervals  much  larger  than  the  damping  band- 
width of  the  system. 
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kx 

/s. 

w 2 

The  magnitude  of  I-;;-!  is  indicated  by  the  shading. 

q 

The  maxima  of  this  function  follow  the  dispersion 
curves  of  the  coupled  system  modes.  Within  each 
interval,  [k  ,k  +dk  ],  the  "spectral  mass"  is  first 

XXX 

resolved  into  the  portions  attributable  to  each 
mode  and  is  then  assigned  to  the  corresponding  dis- 
persion curve.  Finally  each  mass  is  projected  onto 
the  f line  so  as  to  be  distributed  within  the 
intervals  df^  and  df2* 


FIG.  2.11  ILLUSTRATION  OF  THE  APPROXIMATION  USED 
IN  CALCULATING  THE  INTEGRATED  SPECTRUM 
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^ e 


When  the  above  procedure  is  carried  out,  it  is  found 
that  the  integrated  spectrum  can  be  expressed  in  terms  of  the  dis- 
persion relations  of  the  uncoupled  rail  and  slab  modes.  The 
full  result  for  light  damping  is  as  follows: 


np(oj)  = z np  (0)) 

n=l  n 


(2.29a) 
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or 


2 2 2 2 

Here  co  /oo  = A and  co  /oo  = A are  the  dispersion  curves  of 
' or  r or  s ^ 

the  uncoupled  rail  and  slab  modes  and  U is  a unit  step  function; 


U(x)  = 


1 , X > 0 
0 , X < 0 


(2.32) 
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The  value  of  i]>  in  Eqs . 2.29  and  2.30  is  to  be  at  y — y > where 

o 

y is  the  distance  of  the  rail  from  the  edge  of  the  slab,  see 
o 

Fig . 2.1. 


Each  term  in  Eq.  2.29c  shows  a peak  at  the  rail  and 

dk 

slab  resonant  frequencies  due  to  the  dependence  on  This 

result  reflects  the  fact  obtained  in  section  2.4,  that  overall 
slab  response  is  maximum  near  rail  and  slab  resonance.  Below 
these  frequencies  each  term  has  a cutoff  due  to  the  fact  that 
bending  waves  do  not  propagate.  Also,  as  might  be  expected  the 
overall  levels  of  II  are  inversely  proportional  to  the  damping. 

Figs.  2.12  and  2.13  show  the  contribution  of  the  first 

mode,  n averaged  over  5 Hz  bands  for  two  different  cases.  The 
^1 

rail  used  in  the  first  case  is  that  shown  in  Fig.  2.3  with 
(jj  /oj  =10.  Here  no  coincidence  occurs  and  the  reduced  spec- 
trum  falls  off  rapidly  for  frequencies  above  the  rail  resonance. 
In  the  second  case,  the  bending  rigidity  of  the  slab  has  been 
increased  by  ^40%  and  a coincidence  occurs  at  '^>200  Hz.  Although 
the  power  levels  of  the  rail  and  slab  peaks  are  not  appreciably 
affected,  above  u)  = the  reduced  spectrum  shows  a wide  hump 

centered  at  ^180  Hz.  The  level  here  is  some  10  dB  above  v/hat  it 
was  in  the  non-coincident  case.  The  contributions  of  the  higher 
modes  show  similar  characteristics. 

From  the  above  we  conclude  that  if  the  excitation 
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NORMALIZED  INTEGRATED  SPECTRUM 
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NORMALIZED  INTEGRATED  SPECTRUM 


FIG.  2.13  INTEGRATED  SPECTRUM  OF  TUNNEL  FLOOR  FORCE  - WITH  COINCIDENCE 
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spectrum  is  broad  band  and  constant  for  all  frequencies,  the 
overall  magnitude  of  the  tunnel  floor  force  is  greatest  near 
the  uncoupled  rail  slab  resonances.  The  occurrence  of  coinci- 
dence can  significantly  deteriorate  slab  effectiveness. 


wheels  on  track  is  not  constant  for  all  frequencies.  Thus,  in 
section  2.7  the  integrated  spectrum  will  be  calculated  for  an 
excitation  spectra  obtained  from  the  micro-roughness  model  of 
the  wheel-rail  interaction  [5] . 

2,6.2  High  Damping 

For  high  values  of  damping  the  response  of  the  constant 
oscillators  is  predominantly  nonresonant.  Thus,  the  integrated 
spectrum  given  by  Eq.  2.29c  will  be  in  error,  since  it  includes 
only  the  contribution  from  resonant  oscillator  vibration.  If  we 
include  the  nonresonant  response,  the  integrated  spectrum  should 
be  a sum  of  two  terms . 


The  spectrum  of  excitation  due  to  the  rolling  of 


resonant 


nonresonant 


np(u)  = np(w)  + 


(2.33) 


where  the  first  term  is  given  by  Eqs.  2.29. 
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To  calculate  the  nonresonant  contribution  to  the  inte- 


grated spectrum,  we  refer  back  to  Eqs.  2.24,  2.26  and  2.28  which 
give  the  integrated  spectrum  in  terms  of  the  wavenumber  components 
for  each  slab  mode. 


oscillator  at  wavenumber  k . As  noted  previously,  the  spectral 
density  for  light  damping  exhibits  large  maxima  at  values  of  k^ 
and  0)  near  those  given  by  the  dispersion  curves  for  the  uncoupled 
rail  and  slab  modes.  For  large  damping  the  spectral  density  does 
not  exhibit  large  maxima  so  that  all  values  of  k must  be  con- 
sidered  in  the  integration  of  Eq.  2.34.  To  proceed  we  can  assume 
that  wavenumber  coincidence  effects  are  not  important  with  regard 
to  the  nonresonant  response  since  they  occur  only  at  a single 
value  of  wavenumber.  In  this  case  we  can  use  a chaining  approach, 
i.e.,  we  can  calculate  the  rail  response  assuming  the  slab  to  be 
rigid  and  then  use  the  computed  rail  response  to  calculate  the 
forces  on  the  slab,  the  slab  response,  and  the  forces  transmitted 
to  the  tunnel  floor.  Following  this  approach  we  find 


00 


(2.34) 


where  (k^  to)  is  the  spectral  density  of  the  response  of  the  n 
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where  is  defined  in  Fig.  2.1.  The  suinination  and  integration 

required  in  Eq.  2.34  can  only  be  carried  out  numerically.  This 

would  represent  a reasonable  approach  in  detailed  studies  of  one 

specific  design.  However,  our  objectives  for  the  present  study 

are  to  develop  engineering  design  guidelines.  Therefore,  we 

proceed  to  simplify  Eq.  2.35  so  that  the  integration  can  be 

carried  out  analytically  in  closed  form.  Toward  this  end  we 

deal  only  with  cases  where  o)  <<  co  and  where  frequencies  are 
^ os  or  ^ 

in  the  range  co  < co  < co 

^ os  or. 


The  rigid  body  slab  resonance,  co^^,  is  typically  below 

125  radians  per  second  (20  Hz)  and  the  rigid  body  rail  resonance, 

CO  , is  above  785  radians  per  second  (125  Hz).  Thus,  the  fre- 
er 

quency  range  to  co^^  encompasses  frequencies  of  greatest  im- 
portance in  calculating  slab  effectiveness,  since  the  vibration 
transmitted  to  nearby  buildings  is  greatest  in  this  frequency 
range.  Furthermore,  we  deal  only  with  cases  where  the  cutoff 
frequency  of  the  n=2  crossmode,  ^2s'  well  above  co^^,-  In  this 
case  we  can  limit  the  summation  in  Eq.  2.34  to  a single  term 
\/ith  n=l  and;  from  Eq.  2.4,  Eg*  2.34  becomes 


55 


~-r  -r  = (-^)^  (1+6  (1+6  L 

S (oo)  m ' s'  ' r 


(2.36) 


2tt 


dk 


2 

r u)  k 2 

(1+  — ) ^ + 
L ^1 

6^ 

r 

(1 

^ oj  2 

w 

6^ 

s 

2 2 

where  we  have  set  = 1.  As  previously  discussed  we  take 

the  damping  of  the  slab,  6^,  to  be  sufficiently  high  that  the 

contribution  to  the  integral  from  values  of  wavenumber  greater 
2 2 hi 

than  (^2*^  '^^or'  negligible.  Thus,  as  long  as  the  condition 
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is  held  and  6^  < 1,  the  integral  can  be  evaluated  as 
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We  see  from  Eq.  2.39  that  the  resonance  frequency  of  the  slab 
should  be  kept  as  low  as  possible  to  prevent  the  transmission 
of  forces  to  the  tunnel  floor  by  nonresonant  slab  motion.  In 
addition,  the  damping  should  also  be  kept  low.  However,  if  it 
is  too  low  the  resonant  slab  response  will  degrade  the  slab 
perfonnance . 

2.7  Integrated  Spectrum  for  Mi croroughness  Excitation 


of  the  rail  by  the  wheels.  Actual  wheels  and  rails  deviate  from 
perfect  smoothness.  Assuming  uninterrupted  rolling  contact,  the 
rails  will  be  excited  by  the  relative  velocity  of  the  contact 
point  of  the  center  of  the  wheel  and  a rail  reference  axis.  The 
resulting  spectrum  of  the  force  on  the  rail  is  [5] 


v/here  z and  z are  the  rail  and  wheel  impedances  respectively, 
r w 

S is  the  train  speed  and  $^(k)  is  the  wavenumber  spectrum  of 
tlie  relative  velocity  of  wheel  and  rail  due  to  their  combined 
roughness.  has  already  been  obtained.  For  frequencies  above 

roughly  50  Hz,  a reasonable  approximation  for  z^  is  found  to  be 


Here  we  consider  a more  precise  model  of  the  excitation 
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where  M is  the  mass  of  the  wheel,  bearing  and  an  equivalent  axle 
mass . 

To  illustrate  the  effect  of  rail  impedance  on  the 
characteristics  of  the  reduced  spectrum  of  tunnel  floor  force, 
we  shall  assume  M to  be  so  large  that  for  all  frequencies  of 
interest, 

S (m)  03^  |z  I ^ $ (k  = o)/S)  (2.42) 

q ^ X ' ^ 

Thus,  the  wheel/rail  interaction  is  modeled  as  a velocity  source. 
The  effectiveness  of  slab  performance  might  now  be  judged  by 
referring  to  4>^(k);  i.e.,  by  evaluating  IIp/$^(k).  The  varia- 
tion  of  I I , as  determined  by  Eq.  2.18,  radically  alters  the 
character  of  the  integrated  spectrum  associated  with  the  reson9.nt 

slab  vibration.  Figs.  2.14  and  2.15  show  the  first  slab  mode 

contribution  to  the  reduced  spectrum  divided  by  <|)^(k)  for  the 

cases  considered  in  Figs.  2.12  and  2.13.  The  presence  of  resonant 

minima  in  lz^|  at  and  completely  eliminates  the  sharp 

peaks  in  n„  due  to  a flat  excitation  spectrum.  The  effect  of 
^1 

coincidence,  however,  remains  the  same,  i.e.,  to  increase  the  re- 
duced spectrum  about  10  dB  over  a broad  frequency  band  centered 
at  ^180  Hz.  Above  that  frequency,  both  spectra  have  a very  wide 
hump  (not  shown  in  the  figures)  with  a maximum  at  '^^500  Hz. 
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NORMALIZED  INTEGRATED  SPECTRUM  LOG 


FIG.  2.14  INTEGRATED  SPECTRUM  FOR  THE  MICROROUGHNESS 

MODEL -CASE  1 
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1 


NORMALIZED  INTEGRATED  SPECTRUM  ( 


FI6.2.15  INTEGRATED  SPECTRUM  FOR  THE  MICROROUGHNESS 
MODEL  OF  WHEEL-RAIL  INTERACTION-CASE  2 
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It  can  be  concluded  that,  in  general,  for  a wheel/ 
rail  roughness  excitation  the  reduced  spectrum  of  the  tunnel 
floor  force  exhibits  no  sharp  resonant  peaks  near  the  rail 
resonant  frequency  or  the  slab  mode  cutoff  frequencies.  As  in 
the  case  of  a flat  excitation  spectrum,  the  occurrence  of  coin- 
cidence still  has  a significant  effect  on  slab  performance. 

2.8  Total  Acoustic  Power  Radiated  by  the  Slab 


Now  consider  the  total  acoustic  power  radiated  by  the 
slab.  Only  those  wavenumbers  smaller  than  the  acoustic  wave- 
number  (k=Go/c)  may  radiate  sound  efficiently.  But  the  signifi- 
cant response  of  rail  and  slab  (looking  at  the  plane)  lies 

a’^ove  the  acoustic  dispersion  curve.  Then,  for  forcing  at  a 
given  o),  the  total  radiated  power  nay  be  approximated  by: 
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where  p and  c are  the  density  and  velocity  of  sound  in  air,  a 

•^o  o rad 
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is  the  radiation  efficiency  of  the  slab,  and  use  has  been  made 
, of  the  modal  expansion,  Eq . 2.9,  and  the  orthogonality  properties 
of  the  slab  modes.  For  a stationary  random  excitation,  the  spec- 
trum of  the  total  power  radiated  is  therefore: 


”rad‘“>  ' ^ , «rad,n<“> 

n=l 


(2.44a) 
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or,  using  Parseval's  theorem: 
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Comparing  this  expression  with  Eq.  2.28  shows  that 


p c o j „ 

w . (to)  = - 0)2  n 

K?  (It62) 


(2.46) 


where  n (o))  is  the  integrated  spectrum  of  the  force  on  the  tunnel 

r 

n 

floor . 


Apart  from  a factor  of  to  the  total  acoustic  power 

radiated  is  proportional  to  the  reduced  spectrum  of  the  force  on 
the  tunnel  floor  calculated  above.  Thus,  the  general  characteris- 
tics of  the  total  acoustic  power  are  the  same  as  those  derived 
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for  the  reduced  spectrum,  except  for  a magnification  by  a factor 
2 

of  m arad.  The  radiation  efficiency,  a , , can  be  taken  to  be 

rad 

equal  to  one  for  frequencies  above  70  Hz  [6] . 


2.9  Comparison  of  Prediction  Procedure  with  State-of- the-Art  Method 

The  analytical  model  developed  in  the  preceding  sections 
allows  us  to  predict  the  vibration  levels  of  the  rail  and  slab  and 
the  dynamic  forces  on  the  tunnel  floor.  It  is  important  now  that 
we  understand  the  significant  differences  between  our  results  and 
those  obtained  using  the  state-of-the-art  procedure. 

The  state-of-the-art  procedure  is  to  model  the  slab- 
track  by  an  equivalent  single-degree-of-f reedom  system  as  shown 
in  Fig.  1.1.  In  this  simple  model  the  mass  is  taken  to  be  the 
mass  per  unit  length  of  the  floating  slab.  It  is  implicitly 
assumed  that  the  mass  per  unit  length  of  the  rail  is  small  com- 
pared to  that  of  the  slab  and  therefore  can  be  neglected.  The 
spring  constant  in  the  single-degree-of-f reedom  model  is  set  to 
be  equal  to  the  slab  supports  and  the  entrained  air  under  the 
slab  per  unit  length  of  track.  With  these  definitions  the  reso- 
nance frequency  of  the  single-degree-of-f reedom  system  is  simply 
the  rigid  body  resonance  of  the  slab,  as  defined  in  Eq.  2.5b. 
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j?inally,  the  damping  coefficient  for  the  single-degree-of-freedom 
model  is  set  equal  to  the  effective  damping  coefficient  per  unit 


length  of  the  slab  track,  6^. 


The  performance  of  the  floating  slab  is  predicted  in 
terms  of  the  force  transmissibility , Tp(o)),  which  is  defined  as 
the  ratio  of  the  amplitudes  of  the  force  applied  through  the 
spring  and  damping  elements  to  the  rigid  support  and  the  force 
applied  to  the  mass  [7].  This  transmissibility,  T„(m),  is  given 
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w..iere  o)  is  the  frequency  in  radians  per  second.  The  damping  co- 
efficient, 6^,  is  a solid  type  of  damping  [8].  For  a viscous 
damping  mechanism  it  is  more  common  to  use  a viscous  damping  co- 
efficient, C , or  a fraction  of  critical  damping  coefficient,  r 
s s 

These  parameters  are  related  by 
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Thus,  if  the  damping  is  viscous  Eq . 2.47  becomes 
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The  transmissibility  for  constant  values  of  C are  shown  in 

s 

Fig.  1.2 


At  high  frequencies,  Eq.  2.47  can  be  written  approxi- 
mately as 

2 2 I' 

T„(u))  ==  (-^)  (1  + 6“)  2 for  0)  >>  (JO  ^ (2.50) 

r cu  S 03 

or  in  the  case  of  viscous  damping  Eq.  2.49  can  be  written 

(jO 

^ for  o » (2.51) 

The  force  transmissibility  for  the  single-degree-of- 

freedom  system  shows  the  following  characteristics . At  low 

frequencies  below  the  rigid  body  resonance  frequency  of  the 

slab,  to  , the  transmissibility  is  one — forces  are  transmitted 
os 

without  reduction  in  amplitude.  At  and  near  the  resonance  fre- 
quency the  transmissibility  is  greater  than  one  so  that  the  forces 
are  increased  in  amplitude  by  the  addition  of  the  slab.  Finally, 
at  high  frequencies  the  force  transmissibility  is  much  less  than 
one  and  decreases  rapidly  with  increasing  frequency.  For  solid 
type  damping  Eq.  2.50  shows  that  the  transmissibility  decreases 
with  the  second  power  of  frequency  and  does  not  depend  on  the  value 
of  the  damping  unless  6^  is  unusually  large;  i.e.,  near  one.  For 
a viscous  damping  mechanism  Eq.  2.51  shows  that  the  transmissibility 
decreases  with  the  first  power  of  frequency.  In  this  case  the 
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transmissibility  also  increases  with  increasing  values  of 
damping.  Thus,  a tradeoff  occurs.  Small  values  of  viscous 
damping  decrease  the  force  transmissibility  at  high  frequencies, 
but  increases  the  amplitude  of  the  peak  in  the  transmissibility 
at  the  slab  resonance  frequency. 

In  the  more  precise  analytical  model  developed  in  this 
report  the  force  transmissibility  must  be  defined  for  each  cross- 
mode and  as  a function  of  both  frequency,  o),  and  wavenumber,  k^; 
see  Eq.  2.25.  Combining  Eqs.  2.24,  2.25,  and  2.35  with  the  re- 
lationship , 

S (k  oj)  = jw  (k  m)  1^  (2.52a) 

w^  X , ' n X , ' 


and 


S (w)  = |q(w) 


(2.52b) 


we  can  find  the  transmissibility  for  each  cross-mode.  In  the 
special  case  of  the  first  mode  with  k^=0  we  find  the  transmissibility 
to  be 
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For  frequencies  below  the  rail  resonance  frequency  and  for 
moderate  to  small  values  of  rail  damping  the  transmissibility 
given  by  Eq.  2.53  agrees  with  that  given  by  the  single-degree- 
of-freedom  model. 

The  transmissibility  for  the  first  mode  and  for  zero 
wavenumber  indicates  the  transmission  of  forces  that  have  a long 
wavelength.  It  is  reasonable  to  argue  that  these  long  wave- 
length forces  are  more  effective  in  exciting  the  tunnel  into 
vibration  than  the  forces  with  short  wavelength,  since  the  tun- 
nel floor  has  a very  high  bending  rigidity.  Based  on  this  argu- 
ment it  is  also  reasonable  to  conclude  that  the  force  transmis- 
sibility given  by  the  single-degree-of-f reedom  model  may  give  a 
good  indication  of  the  transmission  of  vibration  to  the  tunnel. 
Further  study  of  the  tunnel  response  to  different  types  of  force 
excitation  is  needed  to  resolve  the  general  validity  of  this  con- 
clusion . 


The  transmissibility  given  by  the  single-degree-of- 
freedom  model  also  relates  to  the  transmissibility  for  the  inte- 
grated spectrum  of  forces  due  to  nonresonant  slab  vibration. 
Comparison  of  Eqs.  2.39  and  2.50  shows  that  for  high  frequencies 
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where  H (w)  is  the  integrated  spectrum  from  Eq.  2.39  and  T (to) 

F r 

is  the  force  transmissibility  given  by  Eq.  2.50.  For  typical 
conditions,  RE  100  rail  with  a rail  resonance  frequency  of 
200  Hz  (1257  rad/sec),  Eq.  2.54  becomes 

nonresonant  j 

n„(to)  , (1+<S^)  L 

where  L is  the  slab  width  in  feet.  For  single  track  slabs 
with  a width  of  approximately  10  ft.  and  for  typical  values 
of  rail  fastener  damping,  the  integrated  spectrum  divided  by 
the  spectrum  of  the  forces  applied  to  the  rail  at  the  wheel/ 
rail  contact  point  is  very  closely  related  to  the  force  trans- 
missibility given  by  the  single-degree-of-f reedom  model. 

A comparison  of  results  from  the  precise  model  developed 
in  this  report  with  results  from  the  state-of-the-art  single- 
degree-of-f  reedom  model  shows  that  in  a general  case  where  forces 
over  a broad  range  of  wavenumber  contribute  to  the  tunnel  response 
the  simple  single-degree-of-f reedom  model  can  overestimate  the 
slab  performance  (underestimate  the  transmissibility)  by  over 
10  dB.  Thus,  great  caution  must  be  must  be  used  in  applying  the 
results  of  the  simple  model  even  though  its  use  is  valid  under  a 
number  of  conditions. 
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2.10  Extension  of  the  Analysis  to  Floating-Tie  Track 


In  previous  sections  we  have  carried  out  an  analysis 
for  floating  slabs  that  are  infinite  in  the  direction  along  the 
track.  In  reality,  of  course,  the  slab  is  not  infinite  and  in 
some  cases  is  better  modelled  by  a series  of  short  segments. 
Floating  tie  systems  must  be  included  with  these  cases. 

The  floating-tie  track  consists  of  heavy  concrete 
ties  supported  by  resilient  materials  on  the  tunnel  floor.  The 
rail  is  fastened  to  the  tie  by  means  of  conventional  rail 
fasteners.  To  treat  this  system  analytically  we  must  model  the 
ties  by  an  orthotropic  plate  with  the  bending  stiffness  along 
the  direction  of  the  track  being  zero.  Many  of  the  fundamental 
concepts  discussed  earlier  continue  to  be  valid.  To  calculate 
the  cross-mode  resonance  frequencies  we  use  Eqs.  2.4  and  2.5 
with  k^=0.  The  dispersion  curves,  shown  in  Fig.  2.3  for  a float- 
ing slab,  become  straight  lines  with  values  of  frequency  being 
independent  of  wavenumber. 

The  mechanical  oscillator  analogy  shown  in  Fig.  2.6 
can  also  be  used  for  the  floating-tie  track.  Parameters  are 
obtained  by  setting  k^  equal  to  zero  so  that  the  oscillator 
masses  and  stiffnesses  are  no  longer  functions  of  wavenumber. 
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The  assumptions  and  approximations  leading  to  Eqs. 

2.29  and  2.30  are  no  longer  acceptable  for  the  floating-tie 
system.  However,  if  we  restrict  our  attention  to  the  frequency 
range  below  the  rail  resonance  frequency  we  can  use  the  chaining 
approach,  which  leads  to  an  equation  that  is  similar  to  Eq.  2.35. 
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Following  the  same  approach  used  to  derive  Eq.  2.36  we  find  the 
integrated  spectrum  as 
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For  values  of  6 less  than  one  we  find 
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If  we  limit  the  summation  in  Eq.  2.58  to  the  first  mode,  then 
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The  first  term  in  Eq.  2.59  is  identical  to  the  square  of  the 
force  transmissibili ty  obtained  from  the  single-degree-of- freedom 
model.  The  second  term  is  approximately  equal  to  two  for  typical 
floating-tie  systems.  Thus,  we  conclude  that  the  simple  model 
can  be  used  for  floating-tie  systems  as  long  as  the  resonance 
frequencies  of  the  second  and  other  higher  cross  modes  are  above 
the  frequency  range  of  interest.  If  this  is  not  the  case,  then 
Eq . 2.58  can  be  used  to  provide  a measure  of  performance  for  the 
floating-tie  system. 


2.11  Conclusions  for  Analytical  Model 

The  analytical  model  developed  in  this  section  allows 
us  to  calculate  the  ratio  of  dynamic  forces  applied  to  the  tunnel 
floor  to  those  applied  to  the  rail  at  the  wheel/rail  contact 
points.  Ideally,  we  would  like  to  calculate  the  vibration  levels 
in  buildings  near  the  tunnel.  However,  to  do  this  we  require  in 
addition  models  to  calculate  the  dynamic  forces  at  the  wheel/ 
rail  contact  points,  the  vibratory  response  of  the  tunnel,  and 
the  vibration  transmission  from  the  tunnel  to  nearby  buildings. 
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within  the  present  state-of-the-art  these  models  do  not  exist, 
although  some  relevant  work  has  been  carried  out  [2,9]. 

In  spite  of  our  inability  to  calculate  the  actual 
vibration  levels  in  buildings,  we  can  use  the  analytical  model 
developed  here  to  compare  different  floating  slab  designs. 
Designs  which  minimize  the  integrated  spectrum  of  forces  on 
the  tunnel  floor  in  the  important  frequency  range  from  approxi- 
mately 10  to  175  Hz  are  desired. 
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3.  FIELD  STUDY 


3.1  Introduction 

With  the  assistance  of  the  New  York  City  Transit 
Authority,  Cambridge  Collaborative,  Inc.  has  undertaken  a 
detailed  field  study  of  the  floating  slab  track  section  located 
between  the  Parsons  Boulevard  and  the  169th  Street  stations  of 
the  IND  line.  The  purpose  of  the  study  was  to  gather  compre- 
hensive data  that  would  exhibit  the  slab's  performance  in  re- 
ducing the  vibration  transmitted  to  neighboring  buildings. 

Data  collected  consists  primarily  of  vibration  measurements  at 
several  points  on  the  slab,  on  the  tunnel  walls,  and  on  the 
tunnel  floor.  Measurements  of  the  acoustic  levels  in  the  tun- 
nel were  also  taken. 

Because  the  vibration  levels  arising  from  a train 
passage  vary  from  passby  to  passby  due  primarily  to  differences 
in  wheel  condition  and  train  speed,  it  was  necessary  to  make 
simultaneous  recordings  of  the  vibration  signals  at  many  loca- 
tions in  the  tunnel  for  each  of  several  train  passbys.  With 
these  simultaneous  recordings  it  is  possible  to  compare  vibra- 
tion levels  near  the  floating  slab  track  section  with  levels 
near  a conventional  section  of  track  for  the  same  passby.  The 
comparison  leads  to  a direct  evaulation  of  the  floating  slab 
performance . 
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3.2  Tunnel  and  Slab  Construction 


The  tunnel  is  of  a typical  four-track  cut  and  cover 
construction  with  two  express  tracks  at  the  center  and  two 
local  tracks  on  either  side  as  in  Fig.  3.1.  The  floating 
slab  is  located  on  the  inbound  express  track  and  is  composed 
of  four  fifty-foot  long  sections  which  are  connected  together 
by  steel  reinforcing  rods  as  shown  in  Section  AA  of  Fig.  3.2. 

The  slab  is  isolated  from  the  tunnel  structure  by  means  of 
resilient  rubber-encased  glass  fiber  pads  spaced  as  shown  in 
Fig.  3.2.  The  increasing  number  of  pads  at  both  ends  of  the 
slab  provides  a gradual  transition  of  support  from  the  conven- 
tional track  to  the  floating  slab  track.  A cross  section  detail 
of  the  slab  appears  in  Fig.  3.3. 

The  slab  was  constructed  by  removing  the  existing 
finished  floor  slab  (shown  in  Fig.  3.1)  and  laying  down  a welded 
galvanized  steel  pan  form  with  the  resilient  pads  (Consolidated 
Kinetics  Type  KIP-8802)  epoxied  to  the  pan  bottom.  Steel  rein- 
forcing was  installed  and  finally  the  concrete  poured.  Fiber- 
board  spacers  (Consolidated  Kinetics  Perimeter  Isolation  Board) 
were  set  between  the  slab  and  the  tunnel  piers  to  isolate  the 
tunnel  from  any  lateral  motion  of  the  slab.  A soft  mastic  type 
compound  was  then  poured  on  all  edges  and  seams.  The  measured 
resonant  frequency  is  about  16  Hz  which  indicates  an  overall 
resilient  pad  stiffness  of  2.5  x 10^  Ib^/ft/foot  of  slab.  The 
predicted  resonance  frequency  using  laboratory  measurements  of 
pad  stiffness  was  14  Hz  (not  including  the  weight  of  the  train 
and  the  stiffness  of  the  air  enclosed  under  the  slab)  , 


74 


MANHOLE  ACCESS  , STREET 


75 


Figure  3.1  IND  Line  Tunnel  Cross  Section  Showing  Floating  Slab  and  Tunnel  Access. 
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Figure  3.2  Floating-Slab  Location  Isolation  Pad  Arrangement,  and 
Floating-Slab  Section  Interconnection  Scheme. 
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3,3  Measurement  Scheme  and  Purpose 


The  purpose  of  the  measurement  scheme  was  to  determine 
the  floating  slab's  ability  to  reduce  vibration  transmission  to 
neighboring  buildings.  The  concept  of  Insertion  Loss  was  used 
as  a measure  of  the  slab's  performance.  By  Insertion  Loss  it  is 
meant;  the  reduction  in  vibration  at  a given  point  arising  from 
the  addition  of  the  floating  slab  into  the  tunnel-track  system. 

Fig.  3.4  contains  a diagram  of  the  track  layout  and 
measurement  point  locations.  Notice  that  the  measurement  points 
form  two  tunnel  cross  sections:  one  through  the  floating  slab 

track  section  and  the  other  through  a conventional  track  section. 
Measurement  locations  were  selected  such  that  vibration  measured 
at  the  floating  slab  section  would  be  indicative  of  the  vibration 
that  would  be  observed  if  the  floating  slab  were  infinitely  long. 
Similarly,  measurements  at  the  conventional  track  cross  section 
were  taken  sufficiently  far  from  the  floating  slab  that  they 
would  be  indicative  of  measurements  near  an  infinitely  long 
conventionally  designed  track. 

It  should  be  pointed  out,  however,  that  a New  York 
train  is  normally  600  feet  long  and  the  slab  is  only  200  feet 
long,  so  that  at  any  one  instant  in  time  most  of  the  trucks  are 
located  on  the  conventional  track.  As  a result,  the  vibration 
reduction  that  an  infinitely  long  floating  slab  can  provide 
will  not  be  fully  realized  at  buildings  adjacent  to  the  tunnel. 
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One  question  that  arises  is  whether,  in  fact,  our 
measurements  in  the  tunnel  indicate  the  effect  of  an  infinitely 
long  slab.  As  mentioned  above,  most  of  the  trucks  are  on  the 
conventional  track,  however,  it  will  be  subsequently  concluded 
from  an  examination  of  the  data  that  the  slab  is  long  enough 
to  give  the  same  results  as  an  infinitely  long  slab  for  measure- 
ments at  the  center  island. 

Insertion  Loss  is  arrived  at  by  subtracting  the  ac- 
celeration level  measured  at  the  floating  slab  track  section 
from  that  measured  at  the  same  relative  location  of  the  con- 
ventional track  tunnel  cross  section.  The  measure  of  Insertion 
Loss  that  has  been  considered  is  the  difference  in  acceleration 

levels  at  the  center  island,  positions  8 and  5 (Insertion  Loss  = 
fPos  8)  (Pos 

L^  ‘ — L^  ’ ' ) . When  the  vibration  level  at  point  5 is 

3.  cL 

less  than  that  at  8,  the  Insertion  Loss  is  positive,  i.e.,  the 
effect  of  the  floating  slab  is  to  reduce  vibration  transmitted 
to  the  tunnel  floor  and  .to  adjacent  buildings. 

3.4  Ins  trumen  ta  ti on 

Measurements  were  conducted  by  epoxying  accelerometer 
stud  mounts  to  the  appropriate  positions  in  the  tunnel,  mounting 
the  accelerometers,  and  connecting  the  short  microdot  cables  to 
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preamplifiers  which  acted  as  signal  conditioners.  The  amplifying, 
monitoring,  and  recording  equipment  were  located  in  a commercial 
van  parked  next  to  the  access  manhole  on  the  street.  Single  con- 
ductor coaxial  cables  of  lengths  from  75  feet  to  200  feet  joined 
tunnel  located  transducers  with  the  electronics  in  the  van.  En- 
trance to  the  tunnel  was  gained  through  a manhole  at  the  street 
level  which  opened  into  a side  chamber  as  shown  in  Fig.  3.1. 

Since  low  frequency  information  was  of  prime  importance, 
each  acceleration  signal  was  low  pass  filtered  just  prior  to 
amplification  in  order  to  increase  the  signal-to-noise  ratio  in 
the  lower  frequency  range.  Acceleration  levels  below  125  Hz  were 
of  prime  interest  in  view  of  the  slab's  low  spring-mass  resonance 
(16  Hz) . Low  pass  filters  used  were  single  pole  RC  type  and  had 
a 3 dB  down  point  at  180  Hz.  The  Insertion  Loss  down  to  about 
31.5  Hz  has  been  measured  by  Anthony  Paolillo  [10] , Environment 
Staff  Scientist  for  the  NYCTA,  and  is  in  excellent  agreement 
with  that  measured  here  in  the  frequency  range  31.5  Hz  to  125  Hz. 

Because  of  noise  floor  limitations  in  the  electronics, 
accurate  levels  at  lower  frequencies  could  not  be  determined 
from  Paolillo 's  data.  Particular  attention  has  been  paid  to 
this  problem  in  our  measurements  by  selecting  very  low  frequency 
roll-off  devices.  Typically,  Insertion  Loss  spectra  down  to 
0.63  Hz  have  been  determined. 
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Only  the  most  durable  electronic  system  could  be  used 
because  of  the  severe  working  environment  in  which  measurements 
had  to  be  taken.  Much  effort  was  put  into  collecting  those  com- 
ponents that  would  offer  good  low  frequency  resolution  and  high 
signal-to-noise  ratio.  Typically,  B & K 4344  and  4333  acceler- 
ometers were  attached  to  Ithaco  143  and  125  preamplifiers  which 
provided  a low  impedance  source  to  transmit  data  from  the  tunnel 
to  A.C.  coupled  Ithaco  or  D.C.  coupled  DANA  voltage  amplifiers 
in  the  equipment  van  parked  on  the  street.  The  signals  were  then 
recorded  at  15  ips  on  a Honeywell  5600C,  14-channel,  FM  tape 
recorder  which  has  a maximum  linear  input  of  10  volt  rms  and  a 
D.C.  to  10  KHz  flat  frequency  response.  Amplifier  gains  were 
set  to  their  maximum  linear  voltage  output  by  monitoring  passby 
signals  on  a storage  oscilloscope.  When  all  levels  were  set, 
data  was  then  recorded. 


In  addition  to  acceleration  signals,  acoustic  levels, 
photocell  voltage,  and  an  IRIG  "B"  time  code  were  recorded. 
Acoustic  signals  were  monitored  and  recorded  using  one  inch  GR 
microphones  (with  one  inch  wind  screens)  and  GR  1933  Sound  Level 
Meters.  This  acoustic  information,  though  indicative  of  the  ef- 
fect of  the  slab  on  tunnel  acoustic  levels,  was  only  used  as  a 
reference  since  the  acoustic  performance  is  not  of  great  interest 
at  this  time. 
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The  photocell,  mounted  beside  the  track  in  a position 
such  that  the  car  wheels  broke  its  line  of  sight  with  a reflec- 
tor, provided  an  FM  recordable  DC  voltage  shift  upon  each  wheel 
passage.  This  provided  us  with  both  speed  and  position  informa- 
tion . 


3.5  Representative  Passby  Signals 

The  recorded  data  was  processed  at  Cambridge  Collabora- 
tive's laboratory  facilities.  Of  particular  interest  are  the 
filtered  time  histories  that  were  acquired  by  playing  back  re- 
corded signals  into  four  GR  1933  Sound  Level  Meters  and  plotting 
the  decibel  level  on  a Gulton  Four-Channel  strip  chart  recorder. 
The  photocell  logic  output  was  used  to  drive  an  event  marker  on 
the  strip  chart  recorder  which  was  activated  at  the  passage  of 
each  wheel.  Figs.  3.5  show  a set  of  125  Hz  octave  band  level 
time  histories  and  the  associated  photocell  output.  These  time 
histories  are  useful  in  that  they  give  evidence  of  the  slab's 
performance . 

In  certain  time  histories  the  passage  of  neighboring 
truck  pairs  is  very  evident.  Generally,  the  nearer  the  measure- 
ment position  to  the  track,  the  more  distinguishable  are  the 
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associated  time  history  peaks.  The  comparison  of  time  histories 
for  positions  4 and  8 indicates  this  to  be  very  much  the  case. 

Also,  the  ability  of  the  floating  slab  to  distribute  truck  loading 
and  truck  vibration  is  illustrated  through  a comparison  of  time 
histories  of  channels  5 and  8;  five  being  located  on  the  center 
island  near  the  floating  slab  section  of  track  and  eight  being 
near  the  conventional  track. 

The  effect  of  a finite  length  slab  section  can  be  seen 
by  examining  the  acceleration  time  history  for  position  5.  Know- 
ing that  the  photocell  is  75  feet  from  the  approach  end  of  the 
slab,  it  may  be  observed  that  the  contribution  of  a vibration 
source  at  the  approach  end  to  the  level  at  5 is  about  10  dB  below 
that  measured  with  trucks  nearby  on  the  slab,  i.e.,  the  contribu- 
tion of  off-slab  sources  is  negligible.  Similarly,  the  level  at 
position  8,  75  feet  from  the  approach  end  of  the  slab,  is  negli- 
gibly affected  by  sources  at  or  near  to  the  slab  so  that  it  may 
be  inferred  that  the  position  8 vibration  is  not  influenced  by  the 
slab's  presence.  The  conditions  above  have  been  and  must  be  met 
before  a valid  measurement  of  Insertion  Loss,  in  accordance  with 
the  earlier  stated  definition,  can  be  made.  Similarly,  the  dif- 
ference in  acceleration  levels  at  points  1 and  4 constitutes  a mea- 
sure of  slab  Insertion  Loss  with  regard  to  wall  vibration.  However, 
for  these  points,  the  above-mentioned  condition  that  vibration 
should  be  from  either  slab  or  off-slab  sources,  but  not  both,  is 
only  marginally  complied  with.  This  would  tend  to  result  in  a 
slightly  conservative  estimate  of  Insertion  Loss  because  the 
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measured  transmission  from  the  slab  includes  small  contributions 
from  off-slab  sources. 

.1 

I In  plotting  octave  band  spectra,  steady-state  levels 

j I 

: are  interpolated  from  time  histories  and  used  to  evaluate 

various  aspects  of  slab  performance.  Problems  arise,  however, 

I 

I when  very  low  frequency  analysis  using  time  histories  is  at- 
jl  tempted.  The  fast  averaging  time  of  the  GR  1933  Sound  Level 

j Meter  is  too  short  to  allow  satisfactory  confidence  in  the 

I measurement;  on  the  other  hand,  the  length  of  the  train  passby 

II  is  too  short  to  allow  sufficient  time  for  an  average  level  to 

be  reached  using  the  slow  averaging  sound  level  meter  setting. 
Ideally,  a much  longer  t:.  ain  passby  time  is  needed  for  the 
sound  level  meter  to  reach  a steady-state  level.  The  effect 
of  a longer  train  passby  was  obtained  by  storing  an  eight 
- second  portion  of  predominantly  steady-state  passby  vibration 
levels  in  a B & K transient  recorder.  The  transient  recorder 
digitizes  and  stores  a signal,  and  then  outputs  it  repeatedly 
at  an  adjustable  rate.  The  effect  is  to  compress  the  signal 
time  scale  and  raise  its  frequency  content  proportionally. 

For  instance,  eight  seconds  of  data  can  be  stored  in  the  tran- 
, sient  recorder  and  played  back  in  as  little  time  as  0.32  seconds, 
i which  constitutes  a 25:1  time  compression.  Typically,  data  was 

I recorded  at  15  ips  and  played  back  from  the  tape  recorder  into 

the  B & K transient  recorder  at  60  ips,  and  the  transient  re- 
corder was  adjusted  to  give  a 10:1  time  compression,  thus  ef- 
fectively shifting  the  final  playback  frequency  by  a factor  of 
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40  (16  one-third  octaves)  upwards. 


Using  these  time- compressed  signals,  very  low  frequency 
one-third  octave  band  levels  were  measured  using  a GR  1564  Sound 
and  Vibration  Analyzer  set  to  the  slow  averaging  time.  Slow 
averaged  one-third  octave  band  levels  were  then  read  from  the 
GR  1564  meter.  The  statistical  confidence  limits  here  in  pro- 
cessing the  random  vibration  signals  are  affected  by  both  the 
meter  averaging  time  and  by  the  fact  that  only  an  eight-second 
data  sample  is  used.  A discussion  of  this  will  not  be  attempted 
here,  but  the  resulting  confidence  limits  appear  in  Fig.  3.6. 

The  technique  of  performing  analog  filtering  and 
averaging  on  a long  sample  is  identical  to  the  way  real  time 
analysis  is  performed  by  certain  real  time  analyzers.  Among 
these  are  Spectral  Dynamics  and  Federal  Scientific  which  com- 
press the  time  scale,  sweep  filter  and  average  the  data,  and 
display  the  resulting  frequency  spectra.  This  process  can  be 
continuous,  thus  updating  the  spectra  as  the  time  domain  signal 
changes.  Essentially,  the  sweep  filtering  of  a sample  here  is 
simply  repeated  on  the  same  time  sample  rather  than  updating 
it  by  taking  a new  sample. 

Other  real  time  analyzers  do  not  operate  in  this  man- 
ner. These  input  the  time  domain  signal  into  a bank  of  parallel 
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Figure  3.6  Ninety  Percent  Confidence  Limit  as  a Function 
OF  Frequency  for  Averaging  of  an  Eight  Second  Data  Sample. 
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filters  and  a short  sample  is  taken  from  one  randomly  selected 
filter  band,  digitized,  averaged  with  previous  samples  of  the 
same  band,  and  stored.  This  sampling  occurs  continuously,  and 
an  updated  time  averaged  spectrum  is  generated.  Note  that  this 
is  not  truly  real-time  analysis  because  spectral  levels  are  not 
measured  from  the  same  time  domain  signal  sample.  For  steady- 
state  signals  this  is  usually  permissible,  but  for  short  tran- 
sient signals  erroneous  spectra  can  result. 

3.6  Measured  Insertion  Loss 

Figs.  3.7a,  3.7b  and  3.7c  show  acceleration  levels  in 
third-octave  band.:^  at  measuring  points  5 and  8 for  three  different 
passbys.  Notice  that  at  the  floating  slab's  resonant  frequency 
(16  Hz)  the  level  on  the  tunnel  floor  near  the  slab  track  is 
higher  than  that  at  a conventionally  designed  tunnel  section. 
Hence,  the  slab  has  a detrimental  effect  at  low  frequencies. 

Fig.  3.8  plots  the  corresponding  Insertion  Losses. 

Also  appearing  in  Fig.  3.8  is  the  Insertion  Loss  mea- 
sured by  Paolillo  [10]  . His  levels  extend  up  to  1600  Hz  and  show 
good  agreement  with  our  results  down  to  about  31.5  Hz.  The 
departure  from  our  results  below  this  frequency  is  due  to  a low 
signal-to-noise  ratio  in  his  measurements  at  low  frequencies. 

The  Insertion  Loss  measured  by  Paolillo  is  referenced  to  the 
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• CHANNEL  5 ON  CENTER  ISLAND  AT  FLOATING-SLAB  CROSS  SECTION 
V CHANNEL  8 ON  CENTER  ISLAND  AT  CONVENTIONAL  TRACK  CROSS  SECTION 


Figure  3.7a  Third  Octave  Band  Acceleration  Levels  at  Center 
Island  Points  During17:3A  Train  Passby 
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THIRD-OCTAVE  BAND  ACCELERATION  LEVEL  IN  dB  RE  Ig 


0.5  1.0  2 4 8 16  31.5  63  125  250 

FREQUENCY  IN  HERTZ 
17:40  PASSBY 


• CHANNEL  5 ON  CENTER  ISLAND  AT  FLOATING-SLAB  CROSS  SECTION 
V CHANNEL  8 ON  CENTER  ISLAND  AT  CONVENTIONAL  TRACK  CROSS  SECTION 

Figure  3.7b  Third  Octave  Band  Acceleration  Levels  at  Center 
Island  Points  During  17:A0  Train  Passby 
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0.5  1.0  2 4 8 16  31.5  63  125  250 

FREQUENCY  IN  HERTZ 


17:47  PASSBY 

• CHANNEL  5 ON  CENTER  ISLAND  AT  FLOATING-SLAB  CROSS  SECTION 
V CHANNEL  0 ON  CENTER  ISLAND  AT  CONVENTIONAL  TRACK  CROSS  SECTION 

Figure  3.7c  Third  Octave  Band  Acceleration  Levels  at  Center 
Island  Points  During  17: A7  Train  Passby 
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floor  slab  at  the  center  of  the  track.  It  is  understandable 
that  this  measure  of  Insertion  Loss  is  quite  close  to  that  of 
the  center  island  referenced  Insertion  Loss,  but  how  either 
of  these  relates  to  the  wall  referenced  Insertion  Loss  is  not 
well  understood  within  the  present  state-of-the-art. 

Fig.  3.9  shows  the  difference  between  acceleration 
levels  at  points  5 and  8 for  a train  passage  on  the  outbound 
express  track  (adjacent  to  the  floating  slab  track) . These 
measured  levels  indicate  that  the  addition  of  the  floating 
slab  has  altered  the  response  of  the  tunnel  to  forces  generated 
by  a train  on  conventional  track.  This  result  may  be  due  to  the 
removal  of  the  finished  concrete  track  slab  during  installation 
of  the  floating  slab,  as  shown  in  Fig.  3.1.  The  implication 
here  is  that  the  lightened  tunnel  structure  in  the  area  of  the 
floating  slab  results  in  a generally  higher  vibration  level.  In 
any  case,  we  can  use  the  results  shown  in  Fig.  3.9  to  correct 
the  measured  Insertion  Loss  values  so  that  they  are  applicable 
to  the  case  in  which  a floating  slab  is  added  to  a tunnel  floor 
that  is  identical  at  both  the  floating  slab  track  section  and  at 
the  conventional  track  section.  The  corrected  Insertion  Loss 
values  are  shown  in  Fig.  3.10.  For  this  case  the  floating  slab 
has  little  effect  on  vibration  levels  at  frequencies  below  the 
rigid  body  resonance  frequency  of  the  slab.  At  and  near  the 
resonance  frequency,  the  vibration  levels  are  increased  by  the 
addition  of  the  floating  slab.  And,  finally,  at  high  frequencies 
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(P0S.8)  .(P0S.5) 


0.5  1.0  2 4 8 16  31.5  63  125  250 


FREQUENCY  IN  HERTZ 

□ DIFFERENCE  IN  ACCELERATION  LEVELS 
BETWEEN  POSITIONS  8 AND  5, 

17:34  PASSBY  ON  OUTBOUND  EXPRESS  TRACK 

CHANNELS  ON  CENTER  ISLAND  AT  FLOATING  SLAB  CROSS  SECTION 
CHANNELS  ON  CENTER  ISLAND  AT  CONVENTIONAL  TRACK  CROSS  SECTION 

Figure  3.9  Acceleration  Level  Difference  Between  Points  8 
AND  5 During  The  Passage  of  an  Outbound  Express  Train. 
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VIBRATION  LEVELS  -• ►VIBRATION  LEVELS 

INCREASED  BY  ADDITION  REDUCED  BY  ADDITION 

OF  THE  FLOATING  SLAB  OF  THE  FLOATING  SLAB 


FREQUENCY  IN  HERTZ 


FIG.  3,10  INSERTION  LOSS  CORRECTED  FOR  DIFFERENCES 
IN  TUNNEL  STRUCTURE  AT  THE  TWO  MEASUREMENT  POINTS 
ON  THE  CENTER  ISLAND 
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VIBRATION  LEVELS  — ^VIBRATION  LEVELS 

INCREASED  BY  ADDITION  REDUCED  BY  ADDITION 

OF  THE  FLOATING  SLAB  OF  THE  FLOATING  SLAB 


the  slab  acts  to  significantly  reduce  the  tunnel  floor  vibration 
levels.  This  behavior  is  very  similar  to  that  of  a single-degree- 
of -freedom  system,  see  Fig.  1.1.  A discussion  and  comparison  of 
the  field  results  with  prediction  is  presented  in  section  4. 

The  topic  thus  far  has  been  the  Insertion  Loss  measured 
at  the  center  island.  There  is  some  question  as  to  whether  this 
Insertion  Loss  is  the  same  as  that  which  might  be  measured  at  the 
tunnel  wall.  The  Insertion  Loss  based  on  measured  vibration  levels 
at  two  points  on  the  tunnel  wall — one  near  the  floating-slab  track 
and  one  near  the  conventional  track--is  shown  in  Fig.  3.11.  Be- 
cause of  transducer  failure  we  were  forced  to  use  a transducer  and 
amplifier  at  one  cf  the  tunnel  wall  positions,  channel  4,  that  did 
not  have  adequate  low  frequency  performance.  Thus,  the  Insertion 
Loss  is  plotted  only  for  frequencies  above  20  Hz.  However,  it  is 
important  to  note  that  the  Insertion  Loss  for  wall  vibration  (for 
f > 20  Hz)  is  comparable  to  the  "corrected"  IL  for  center  island 
vibration  (Fig.  3.10).  In  spite  of  the  relatively  short  length  of 
the  slab  it  appears  to  be  effective  in  reducing  vibration  levels 
both  of  the  tunnel  floor  and  the  tunnel  walls. 

The  question  as  to  whether  the  floor  vibration  or  wall 

vibration  is  more  significant  in  determining  the  vibration  of 

nearby  buildings  is  difficult  to  answer.  We  point  out,  however, 

that  the  floor  vibration  levels  are  at  least  20  dB  greater  than 

wall  vibration  levels,  as  shown  in  Fig.  3.12.  This  suggests 

that  the  floor  vibration  levels  may  be  more  important  with  regard 

to  vibration  transmission  to  buildings  than  the  wall  vibration. 

98 


THIRD-OCTAVE  BAND  INSERTION  LOSS  IN  dB 


FREQUENCY  IN  HERTZ 


^ 17--34PASSBY 
□ 17:40  PASSBY 


INSERTION  LOSS  AT 

(P0S.4)  (POS.l) 

WALL(Lo  -Lo  ) 


CHANNEL  4 ON  WALL  AT  CONVENTIONAL  TRACK  CROSS  SECTION 
CHANNEL  1 WALL  POSITION  AT  FLOATING  SLAB  CROSS  SECTION 


FIG.  3.11  INSERTION  LOSS  FOR  VIBRATION  OF  THE  TUNNEL  WALL 
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VIBRATION  LEVELS  ► VIBRATION  LEVELS 

INCREASED  BY  ADDITION  REDUCED  BY  ADDITION 

OF  THE  FLOATING  SLAB  OF  THE  FLOATING  SLAB 
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0.5  1.0  2 4 8 16  31.5  63  125  250 


FREQUENCY  IN  HERTZ 
17:40  PASSBY 


A CHANNEL  1 WALL  POSITION  AT  FLOATING  SLAB  CROSS  SECTION 
• CHANNEL  5 ON  CENTER  ISLAND  AT  FLOATING  SLAB  CROSS  SECTION 

Figure  3.12  Third  Octave  Band  Acceleration  Spectra  of  Wall 
AND  Center  Island  Vibration  (Points  1 and  5). 
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3.7  Vibration  Due  to  Static  Load 


The  measured  vibrational  acceleration  levels  at 
frequencies  below  approximately  10  Hz  are  quite  high.  The 
mechanism  leading  to  this  low  frequency  vibration  is  believed 
to  be  the  passage  of  the  static  train  load  by  the  measurement 
point.  Confirmation  of  this  explanation  was  obtained  by 
examining  the  spectral  content  of  the  photocell  logic  signal. 

This  signal  is  a square  wave  which  fluctuates  between  two  voltages 
based  on  whether  the  light  beam  is  blocked  by  a train  wheel  or 
not.  The  one- third  octave  band  spectrum  of  the  photocell  signal 
is  shown  in  Fig.  3.13.  Also  contained  here  are  the  one-third 
octave  band  acceleration  spectra  for  three  passbys  at  points  7 
and  9.  Notice  that  the  peaks  in  the  acceleration  spectra  are 
clearly  identifiable  with  those  in  the  photocell  signal  spectra. 

By  knowing  train  speed  and  truck  wheel  bases,  specific  passage 
frequencies  can  be  identified.  The  three  that  most  significantly 
stand  out  are  the  car  passing  frequency,  the  passage  frequency 
associated  with  the  neighboring  trucks  on  adjoining  cars,  and 
the  wheel  passage  frequency  associated  with  the  distance  between 
axles  on  a truck.  These  passage  frequencies  are  identified  in 
Fig.  3.13  for  those  trains  involved.  Data  presented  are  from 
trains  which  have  a truck  wheel-to-wheel  distance  of  6'  - 10", 
a car  length  of  60  feet,  and  neighboring  truck  distance  of 
15  feet.  These  dimensions  are  outlined  in  Fig.  3.14. 
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17:34  PASSBY 

Q CHANNEL  9 AT  FLOOR  BENEATH  CENTER  OF  FLOATING  SLAB 
• CHANNEL  7 ON  FLOATING  SLAB  AT  CENTER 
y CHANNEL  10  (PHOTOCELL) 

A WHEEL  PASSAGE  FREQUENCY  = 8.6  Hz 
• NEAR  TRUCK  PAIR  PASSAGE  FREQUENCY  = 3.9  Hz 
A CAR  PASSAGE  FREQUENCY  = 0.98  Hz 

VELOCITY  = 58.8  FT/ SEC 

Figure 3.13a  Third  Octave  Band  Acceleration  Levels  at  Center  of 
Floating  Slab  (7)  and  Beneath  Floating  Slab  (9)  for  17:34  Passby 
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0.5  1.0  2 4 8 16  31.5  63  125  250 

FREQUENCY  IN  HERTZ 

17=40  PASSBY 

• CHANNEL  9 AT  FLOOR  BENEATH  CENTER  OF  FLOATING  SLAB 
^ CHANNEL  7 ON  FLOATING  SLAB  AT  CENTER 

Figure  3.13b  Third  Octave  Band  Acceleration  Levels  at  Center  u: 
Floating  Slab  (7)  and  Beneath  Floating  Slab  (9)  for  17; 40  Passey 


103 


FREQUENCY  IN  HERTZ 
17 :47  PASSBY 

• CHANNEL  9 AT  FLOOR  BENEATH  CENTER  OF  FLOATING  SLAB 
□ CHANNEL  7 ON  FLOATING  SLAB  AT  CENTER 
O CHANNEL  10  (PHOTOCELL) 


^ WHEEL  PASSAGE  FREQUENCY  = 9.35  Hz 
▼ NEAR  TRUCK  PAIR  PASSAGE  FREQUENCY  = 4.25  Hz 
^ CAR  PASSAGE  FREQUENCY  = 1.06  Hz 
VELOCITY  = 63.8  FT/ SEC 

Figure  3.13c  Third  Octave  Band  Acceleration  Levels  at  Center  of 
Floating  Slab  (7)  and  Beneath  Floating  Slab  (9)  for  17:A7  Passby 
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Figure  3.14  Car/Wheel  Arrangement 


It  may  then  be  concluded  that  the  1.0  Hz  peak  in  the 


acceleration  spectra  results  from  the  passage  of  the  static 
train  load.  The  importance  of  this  peak  with  reference  to  propa- 
gation to  neighboring  buildings  can  be  determined  by  examining 
signals  coming  directly  from  the  recorded  acceleration  channels. 
Fig.  3.15  contains  0.5  Hz  to  25  Hz  band  filtered  acceleration 
signal  time  histories.  They  are  predominantly  of  very  low  fre- 
quency and  basically  occur  only  when  the  train  is  in  the  immedi- 
ate vicinity  of  the  accelerometer  positions.  The  implication 
here  is  that  this  low  frequency  excitation  is  non-propagating 
and  is  a decaying  near  field  effect  characteristic  of  rolling 
contact  deflection. 

Since  the  wall  positions  are  some  distance  from  the 
track,  it  is  expected  that  the  1.0  Hz  spectral  peak  would  be 
less  evident.  Comparison  of  wall  and  center  island  acceleration 
spectra  shown  in  Fig.  3.12  bears  this  out. 

It  should  also  be  pointed  out  that  the  measurement  of 
low  passage  frequency  levels  is  not  subject  to  the  confidence 
limits  of  Fig.  3.6.  These  confidence  limits  apply  to  a station- 
ary random  process,  whereas  the  passage  frequency  levels  are 
from  a deterministic  excitation  mechanism;  i.e.,  the  passage  of 
various  specific  train  components,  and  are  exact  within  limits 
of  small  transducer  and  electronics  errors.  Only  the  wheel/rail 
roughness  is  assumed  to  be  random  and  the  confidence  limits  men- 
tioned are  applicable  only  insofar  as  Gaussian  Statistics  validly 
describe  the  excitation. 
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ACCELERATION  IN  g’s  (0.5  Hz -25  Hz  BANDPASS  FILTERED) 


CHANNEL  5 ON  center  island  AT  FLOATING  SLAB  CROSS  SECTION 


CHANNEL  8 ON  center  island  AT  CONVENTIONAL  track  CROSS  SECTION 


2 SEC 


CHANNEL  3 ON  FLOATING  SLAB,  AT  TRACK  Q.  15'  FROM  SECTION  END 


CHANNELS  AT  FLOOR  BENEATH  CENTER  OF  FLOATING  SLAB 


Figure3.15  AC  Acceleration  Signals  Measured  at  Points  5.  8,  3 and  9. 
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3.8  Slab  Support  Damping 


The  damping  of  the  slab  supporting  elements  plays  an 
important  role  in  determining  the  forces  transmitted  to  the  tun- 
nel floor.  Prior  to  construction  of  the  floating  slab,  tests 
were  run  to  determine  the  stiffness  and  damping  of  the  support 
pads  [11]  . The  results  showed  the  pad  stiffness  to  be  in  the 
range  10,740  to  30,650  Ibs/in  depending  on  load.  The  stiffness 
for  load  corresponding  to  the  slab  weight  is  approximately  12,000 
Ibs/in.  The  damping  of  the  pads  was  in  the  range  2 to  2.4%  of 
critical  damping. 

The  meai>urcd  values  of  Insertion  Loss  show  a negative 
peak  of  -6.3  dB  at  the  slab  resonance  frequency.  From  the 
analytical  model  this  would  indicate  a damping  that  is  30%  of 
critical  damping. 

The  discrepancy  between  these  two  damping  values  may 
be  explained  by  the  addition  of  Consolidated  Kinetics  Transi- 
dampers  10  ft.  on  centers  at  the  edges  of  the  slab.  These  were 
installed  after  the  design  shown  in  Figs.  3.2  and  3.3  and  before 
our  measurements.  The  damping  of  the  Transidampers  has  not  been 
measured  as  part  of  this  study  and  data  were  not  available. 
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3.9  Conclusions  for  the  Field  Study 


The  floating  slab  in  NYC  is  quite  effective  in  reducing 
the  forces  transmitted  to  the  tunnel  floor  at  frequencies  above 
25  Hz.  However,  the  slab  has  a slightly  detrimental  effect  near 
its  resonant  frequency  (16  Hz).  Fortunately  this  frequency  is  at 
a point  where  acceleration  levels  of  the  tunnel  are  at  a minimum, 
as  shown  in  Figs.  3.7  and  3.12.  The  Insertion  Loss  for  vibration 
of  the  tunnel  wall  is  similar  to  that  for  vibration  of  the  center 
island  above  20  Hz.  Insertion  Loss  for  wall  vibration  at  lower 
frequencies  was  not  determined  because  of  limitations  in  the  ac- 
celerometer electronics  at  the  wall  location  near  the  conventional 
track  section.  However,  at  these  low  frequencies  the  wall  levels 
at  the  floating  slab  track  section  were  shown  to  be  much  lower 
than  the  floor  acceleration  levels.  The  question  arises  as  to 
which  contributes  more  to  vibration  transmission  to  neighboring 
buildings.  In  any  case  the  Insertion  Loss  measured  at  the  center 
island  should  be  representative  of  the  vibration  isolation  effect 
of  an  infinite  slab. 

The  effect  of  the  floating  slab  becomes  less  evident  as 
the  measurement  position  moves  further  from  the  track.  The  indica- 
tion is  that  the  vibration  isolation  of  the  slab  is  not  fully 
realized  in  neighboring  buildings  because  it  is  shorter  than  the 
trains.  Further  measurements  would  be  necessary  in  order  to  for- 
mulate a slab  length  criteria  for  reducing  vibration  transmission 
to  neighboring  buildings. 
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The  vibration  levels  at  and  near  1 Hz  are  quite  high. 
These  levels  are  believed  to  be  due  to  the  passage  of  the  static 
load  of  the  train  at  the  car  passing  frequency,  which  is  calcu- 
lated to  be  0.98  Hz.  Vibration  levels  on  the  tunnel  wall  at  and 
near  1 Hz  are  approximately  30  dB  below  the  levels  on  the  center 
island  and  are  only  slightly  higher  than  levels  at  higher  frequen- 
cies. The  importance  of  these  very  low  frequency  (1  Hz)  vibrations 
in  terms  of  vibration  transmission  to  neighboring  buildings  has 
not  been  determined. 
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4.  CONCLUSIONS  AND  DESIGN  RECOMMENDATIONS 


The  work  described  in  this  report  has  focused  on  a 
model  to  predict  the  vibration  isolation  (or  reduction)  pro- 
vided by  a floating  slab  or  tie  track.  In  Section  2 we  have 
shown  the  development  of  a model  which  allows  us  to  relate  the 
forces  acting  on  the  tunnel  floor  to  the  forces  applied  at  the 
wheel/rail  contact  points.  Within  the  limitations  of  the  pres- 
ent study,  we  did  not  extend  the  model  so  that  the  response  of 
the  tunnel  floor  could  be  predicted  or  so  that  the  forces  on 
the  tunnel  floor  could  be  predicted  without  knowing  the  forces 
on  the  rail.  Because  of  this  limitation  we  could  not  use  data 
from  the  field  study  to  verify  the  analytical  model.  We  had 
hoped  to  use  rail  and  slab  vibration  data  together  with  data 
on  rail  fastener  and  slab  support  pad  stiffnesses  to  calculate 
the  forces  acting  on  the  slab.  However,  we  were  unable  to  do 
this  because  we  could  not  collect  rail  vibration  data  at  a suf- 
ficient number  of  locations  so  as  to  specify  the  spatial  dis- 
tribution of  the  vibration  levels. 

A basic  question  that  must  be  answered  before  the 
analytical  model  can  be  applied  is  how  the  tunnel  responds  to 
the  dynamic  forces  applied  by  the  slab  supports.  If  we  assume 
the  tunnel  floor  to  be  much  stiffer  than  the  slab  then  the 
dynamic  forces  with  a short  wavelength  in  space  (high  wavenumber) 
will  be  less  effective  in  exciting  the  tunnel  than  the  forces 
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with  long  wavelength.  In  the  limit,  we  can  assume  that  forces 
due  only  to  vibration  of  the  slab  in  its  1st  crossmode  (n=l) 


and  for  k^=0  are  important.  In  this  case  the  analytical  pre- 
diction model  yields  the  result  that  the  force  transmissibility , 
the  ratio  of  the  amplitudes  of  the  force  on  the  tunnel  floor  to 
the  force  on  the  rail,  is  given  approximately  by 


v/here  6^  is  the  slab  support  damping  loss  factor,  f^^  is  the 
resonance  frequency  of  the  rigid  slab  on  the  support  stiffness, 
and  f is  the  frequency  in  Hertz.  This  result  is  identical  to 
that  given  by  the  simple  single-degree-of-f reedom  model  shown 
in  Fig . 1,1. 

The  expression  given  by  Eq,  4.1  is  for  the  case  v/hen 
the  damping  of  the  slab  supports  is  of  the  'solid'  type.  In 
many  cases  the  damping  mechanism  is  better  modeled  by  a viscous 
damping  factor,  which  is  related  to  the  solid  type  damping 

loss  factor  by  [8] 


(1  + 6^) 


(4.1) 


s 


(4.2) 
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The  force  transmissibility  for  viscous  damping  is  found  by 
substituting  Eq.  4.2  in  Eq.  4.1. 

In  the  most  general  case  the  damping  coefficient 

is  frequency  dependent  and  cannot  be  modeled  by  either  the 

solid  type  or  the  viscous  type  over  the  entire  frequency 

range  of  interest.  In  this  case,  we  let  6 or  r be  a 

s s 

function  of  frequency.  The  force  transmissibility  is  still 
obtained  from  Eq.  4.1. 

The  results  of  our  field  measurements  in  New  York 
support  the  validity  of  the  assumption  that  the  force  trans- 
missibility given  by  the  single-degree-of-f reedom  system  can 
be  used  to  predict  the  performance  of  a typical  slab  in 
reducing  the  tunnel  vibration.  However,  it  must  be  pointed 
out  that  in  very  lightly  damped  systems,  the  vibration  in- 
duced by  higher-order  slab  crossmodes  and  at  higher  wave- 
numbers  may  degrade  the  slab  performance.  The  extent  of 
the  degradation  can  be  formed  by  using  the  analysis  presented 
in  Section  2. 

Floating  tie  systems  should  give  the  same  or  better 
performance  as  given  by  long  slabs,  since  wave  propagation 
down  the  slab  is  inhibited.  If  we  consider  only  the  rigid 
body  response  of  the  ties,  i.e. , no  bending  deformation  across 
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the  tracks,  then  the  performance  of  the  float ing-tie  system 
can  also  be  predicted  using  the  results  from  the  single- 
degree-of-f reedom  model. 

A number  of  design  recommendations  can  be  made  based 
on  the  results  of  this  study.  The  stiffness  of  the  floating 
slab  support  elements*  should  be  selected  so  that  the  deflection 
of  the  rail  under  anticipated  train  loads  is  within  limits  set 
by  safety  and  ride  comfort.  The  mass  per  unit  length  of  the 
floating  slab  or  floating-tie  system  should  be  set  so  that  the 
resonance  frequency  of  the  slab  on  ties  on  their  supporting 
elements  is  in  the  range  8 to  25  Hz.  The  performance  of  the 
8 Hz  system  in  reducing  vibration  will  be  5 to  10  dB  better 
than  that  of  the  25  Hz  system. 

The  resonance  frequency  should  not  be  set  below 
8 Hz  because  of  the  possible  amplification  of  the  low  frequen- 
cy vibration  which  is  caused  by  the  passage  of  the  static 
train  load.  On  the  other  hand  the  resonance  frequency  should 
not  be  set  above  25  Hz  because  of  possible  amplification  of 
the  vibration  in  the  important  31.5  Hz  octave  band. 

As  a general  rule  the  mass  per  unit  length  of  the 
floating-tie  system  must  be  greater  than  that  of  a floating 
slab.  This  is  due  to  the  fact  that  the  rail  deflection  under 

*including  the  effective  stiffness  of  the  air  trapped  under  the 
slab 
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the  train  load  is  reduced  by  the  bending  stiffness  of  the 
slab,  thereby  allowing  the  slab  supports  to  have  less  stiff- 
ness. It  follows  that  the  mass  required  to  achieve  a given 
resonance  frequency  will  be  less  for  the  slab  than  for  the 
tie  system. 

Damping  of  the  slab  supports  also  plays  an  important 
role  in  determining  the  vibration  reduction.  The  damping 
should  be  as  high  as  possible  to  limit  the  amplification  of 
the  vibration  at  the  resonance  frequency.  However,  in  many 
viscoelastic  materials  the  damping  modulus  is  frequency  depen- 
dent and  is  best  modeled  by  a viscous  damping  mechanism.  Then, 
as  shown  in  Fig.  1.2,  high  values  of  damping  decrease  the  per- 
formance of  the  floating  slab  or  tie  system  at  frequencies  above 
the  resonance  frequency.  The  designer  should  strive  to  use 
materials  for  which  the  damping  modulus  and  the  Young's  modulus 
do  not  increase  with  increasing  frequency. 
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APPENDIX  A: 


REPRESENTATIVE  VALUES  OF  PARAMETERS 


Rail : 

= 34  Ib/ft 

Ej,  = 4.32  X 10^  Ib/ft^  = 30  X 10^  Ib/in^ 

= 49  in'*  = 2.363  x 10“^  ft'* 

Slab; 

h = 1 ft 

L = 10  ft 

Pg  = i.O  Ib/ft^ 

Eg  = 3.5  X 10^  Ib/in^  = 5.04  x 10®  Ib/ft^ 

(jo  = 62.8  rad/sec. 

os 


In  addition  the  values  of  various  subsidiary  parameters  are  a 
follows : 


K 


^ = .0410  ft  '* 


1 “ E I 
r r 

2 

P w /I 

<2  = = . 0430  ft  '* 


E I = 1.02  X 10^  Ib-ft^ 
r r 


D = 4.3  X 10  Ib-ft 
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. 0226 


e = = 

A 

T = 2/Kt  k„3  = .0217  ft 
1 2 
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APPENDIX  B:  ORTHOGONALITY  OF  THE  SLAB  MODES 


The  free,  undamped  motion  of  the  slab  is  governed  by 


D[^  + 2 

9x 


2 2 
9x  9y  9y 


+ w + kgW  + pg  = 0 


(B.l) 


9t 


The  free  end  conditions  on  the  edges  of  the  slab  imply  the 
following  boundary  conditions 


2 2 3 3 

9w,  9w  ,,  9w,,^^9w  ,, 

+ y o = 0,  y + (2-y)  5 ^ = 0;  y=0,L  (B.2, 


9y 


9x 


9y 


2 2 
9x  9y 


When  Eq.  2 3 is  substituted  into  *-be  above,  it  is 
found  that  the  ip  must  satisfy 


4 

d ip 
dy 


==  dy2 


4, 

Y ^ 
n n 


(B,3' 


with  the  boundary  conditions: 


d^ip 


n 


dy' 


9 

yk  \p  , ^ 

X n , 3 


= (2-y)  k 


dy' 

for  y = 0 , L 


2 di|; 
X dy 


(B.4) 


where 


4 4,4 

Y = a - k 
' n n X 


(B.5) 
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and  a is  given  in  terms  of  the  frequency  by  Eq.  2.5a. 

The  essential  question  concerns  the  orthogonality  of 

(Y/k  ) for  a given  k . That  is , k is  a parameter 
which  is  fixed  during  the  course  of  this  discussion.  The  eigen- 
value problem  represented  by  Eq.  B.3,  and  Eq.  B.4  thus  determines 

4 4 

discrete  values  of  y m terms  of  k , one  y for  each  distinct  ip  . 

' n X ' n ^n 

To  investigate  orthogonality,  consider  the  quantity 


d"^ip 


dy 


d^ip 


n 


dy' 


m 


dy 


where  n m. 


Integrating  by  parts,  one  has 


L 

f 

I 

0 


n 


dy 


d^ip 

^^dy  = 
dy 


2k 


X 


2 

0 


- [ij;'  (i|;'-2k  : 

m n x^n 


n m 


[ll; 

n m 


L 

0 


ilj  (\h'" 


2k  ip")  dy 
x^m 


(B.6) 


where  the  primes  denote  differentiation  with  respect  to  y.  Sub- 
stituting the  boundary  conditions  (Eq.  B.4)  for  into  Eq.  B.6 


gives : 


121 


(B.7) 


x^n  ^ m x^m 


dy 


= [ ip'  (lb"  -uk^ilJ  (iJj"' 


k^  (y-2)i|;'  ) ]!; 

X 0 


First,  note  that  from  Eq.  B.3  that  the  left  hand  side  becomes 


(Y 


4 

n 


Note  also  that  ^ must 
m 

so  that  the  right  hand 
Hence,  we  obtain 


n 


m 


satisfy 
side  of 


dy 


the  boundary  conditions  Eq.  B.4 
Eq.  B.7  is  identically  zero. 


- ^m) 


L 

f 

^ 0 


dy  = 0 


(B.8) 


which  states  that  for  each  value  of  , the  corresponding  set 
of  eigenfunctions  of  Eq.  B.3  are  orthogonal. 
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APPENDIX  C:  PARTIAL  DIFFERENTIAL  EQUATIONS  OF  MOTION  AND  THEIR 

SOLUTION 

Here  we  set  forth  the  mathematical  background  for  the 
work  presented  in  Section  2.  The  basic  idealizations  used  for 
this  calculation  are  given  on  Page  14.  In  addition  we  limit 
our  consideration  here  to  the  case  in  which  the  rail  is  excited 
by  a point  force,  q(t).  With  these  assumptions  and  the  coordi- 
nates and  conventions  shown  in  Fig.  2.1,  the  equations  of  motion 
for  the  rail  displacement,  ?,  and  slab  displacement,  w,  are 

4 2 

E I —4  + T + K*  [?  - w(x,y  ,t)]  = q(t)  6 (x)  (C.la) 

^ ^9x^  ^ 8t^  ^ ° 

_ 2 2 2 

w + K*w  + pg  = K*  [(?  - w)]  6 (y  - y^)  (C.lb) 

i 3x  9y  : 

!_  j 

The  stress-free  condition  on  the  edges  y = 0,  L of 
the  slab  give  rise  to  the  following  boundary  conditions: 

2 2 

^ f =0  at  y = 0,L  (C.2a) 

3y^  3x 


^ + (2  - y)  =0  at  y = 0,L  (C.2b) 

3y  3x  3y 

where  y is  Poisson's  ratio  for  the  slab  material. 
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The  total  system  response  is  formulated  by  expressing 


the  equations  of  motion,  Eqs.  C.l,  in  terms  of  the  Fourier 
transforms  with  respect  to  x of  the  rail  and  slab  displacement. 


? (k 

X 


I -jk  X 

t)  = } dx  e ^ 5 (x,  t) 

j 

— 00 


(C.3a) 


; -]k  X 

w(k^^y,t)  = dx  e ^ w(x,y,t) 


(C.3b) 


Then,  we  note  that  w may  be  generally  expressed  as 


w(k  y,t) 


CO 


n=l 


w 

n 


(k^^t)  *^(y) 


(C.4) 


where  w are  the  modal  coordinates  of  the  slab  motion.  Sub- 
n 

stituting  Eq.  C.4  into  Eqs.  C.l,  multiplying  both  sides  of 
Eq.  C.lb  by  (for  some  m)  and  integrating  over  y = C,L 

gives  us  the  equations  of  motion  for  the  rail  displacement 
and  the  slab  modes.  Finally,  neglecting  coupling  terms  due 
to  damping , we  find 


K ^ 00 

^ 5 +[k^  + <^(l-j6^)l  5-  E 

0)  m=l 

or 


(C . 5a) 


w + [k^  + ^ K„(l-j6  ) + 3K^o(l-j<S  ) 4^  (y  ) ] w 

2 n X 'n  2 2 s 2 r ^n  -^o  n 

0)  0) 
or  or 


^2®  *n‘yo>  ^ 

m?^n 


(C. 5b) 
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where  Q(t)  = q(t)/EI 


and  3 = 


LP. 


For  floated  slab  systems  it  is  realistic  to  assume  both  and 
are  much  less  than  1 and  that  3 is  much  less  than  1.  For  3 
small,  the  total  system  behavior  is  determined  primarily  by  the 
relative  positions  of  uncoupled  rail  and  slab  dispersion  curves 
shown  for  a typical  case  in  Fig.  2.3. 


Explicit  solutions  to  Eqs.  C.5  may  be  found  by  ex- 
pressing these  equations  in  terms  of  the  space-time  transforms 
of  C and  w.  We  define 


oo 

C ='  dt  e^""^  C(k  t) 

, ^ 9 

— CO 

oo 

w = : dt  e^^^  w (k  t) 

n ■ n X , 

00 

oo 

Q = [ dt  Q(t) 

i 

oo 

In  terms  of  these  quantities,  Eqs.  C.5  become 


(C. 6a) 


(C.6b) 


(C.6c) 
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00 


(C.7a) 


A? 


K S 
^ m=l 


w 

m 


b w 
n n 


^2®’"n‘yo'  ^ 

m=l 


(C.7b) 


where 


j 6 K 

r 


1 


b 

n 


ji<o  C 
2 sn 


os 


or 


sn 


6 

s 


+ 63 

r 


2 

to 


or 


os 


'P 


2 

n 


(C. 8a) 


(C. 8b) 


(C. 8c) 


This  infinite  set  of  equations  may  be  solved  explicitly 

for  and  each  of  the  w . From  these  results  it  can  be  shown 

m 

that  in  the  wavenumber  and  frequency  region  in  which  the  un- 

th 

coupled  rail  dispersion  curve  passes  closer  to  the  n^  slab  mode 

^ /s 

dispersion  curve  than  to  any  other  mode,  then  C and  w^  can  be 

approximately  obtained  by  considering  the  interaction  of  the 
th 

rail  with  the  n mode  alone.  This  mode  by  mode  interaction  ap- 
proximation allows  relatively  simple  solutions  to  be  obtained. 


To  find  the  spatial  characteristics  of  response,  it 
is  assumed  that  only  the  first  slab  mode  is  significantly 
excited  although  qualitative  results  will  hold  in  the  general  case. 
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/\ 

First,  we  neglect  all  quantities  except  ? and  in 


Eqs.  C.7  to  obtain 

/\  A. 

5 = ^ 

(C.9a) 

A A ^ 

where 

(C.9b) 

H = ®i-r 

(C.  10a) 

^W1  “ 

(C.lOb) 

and  where 

^ 2 
V = B^A  - >^2.'^2^^1 

(C.lOc) 

4 4 

A = k - k - j6  Kt 
X r r 1 

(C.lOd) 

2 

a M CA3  ry 

Bi  - kx  - a - :<2  2 

03 

or 

(C. lOe) 

4 

Expanding  the  denominator,  V , in  k , it  is  found  to 
be  factorable  in  the  form; 
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+ X - j A ) (k 
r r X 


4 2 

a + 


- X 


- jA^) 


(C.ll) 


where 


A = 


6 K,  - 

r 1 I 


A os  ^ 

^2  2 
0) 

or 


A^  = ^[A  - + 4k^K23<^^] 


CO 


A = - -=2  «sll 


AA~+~Ak 2?>Ai  “ A 


2 2 ^sl' 

(0 

or 


v4^ 


4 4 2 

A = k - a + 

r 2^1 


Now  consider  the  behavior  of  the  time  transforms  of 
^ and  w^;  i.e. , 


2it 


jk  X ^ ~ 

jdk^  e ^ C, 


2tt 


jk  X 

dkx  e ^ w^ 


(C.12) 


For  a given  (o,  Eqs.  C.IO  and  C.ll  show  that  the  magnitudes  of 
the  integrands  are  sharply  peaked  at  the  intersections  of  un- 
coupled rail  and  slab  dispersion  curves  with  co  = constant. 

The  parameters  A^  and  A^  represent  the  "spread"  of  these  peaks 
about  their  maxima,  while  |A-  2A^|  is  the  separation  be'  =»' 
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If  damping  is 


maxima  and  is  always  greater  than  2/ic^^K^ 
sufficiently  small  that  2 /k  ^ r ^ 3 >>  or  (i.e., 

forces  are  much  larger  than  damping  forces) ; then  the 
forms  of  rail  and  slab  response  are  approximately: 


; A - A 
C _ r 


A 


I - 


A - 2A  r . s 

^ r A - 2A 

Q r 


where 


3k  4) 

— = ^ ^ - I } 

- A - 2A  ^^r  s^ 

Q r 


r 


4lk, 


/r  'Ti 1^1 

e [cos  X + sin  x]  ; 


/2 

~4  I ~4  I 

< 0 & |kj|  >>  A^ 


/2 


I 4/2A 
1 r 


3/4  [(1  + )^  + j(i  - i )^] 

3/4  J2  /T 


'r  ] 


— L_  [j(i  + Ji-)^  - (1  - — )^] 


/T 


/2 


/2 


+ j e 


-X  — [j(i  - — )^  + (1  + — )^]  ; 

/2  ^ /2  /2 


1 L 

4 iT. 


-X  k 
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(j|k, 


4|J  13 


-e 


+ D e 


) ; 


> 0 |k  r >>  A 


i. 


coupling 
time  trans- 

(C.  13a) 
(C.13b) 


(C. 14a) 


(C.14b) 


(C. 14c) 
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and  where 


I = I as  k approaches  a and  A approaches  A (C.15) 

S 3T  3T  JT  S 


Both  rail  and  slab  motions  are  seen  to  consist  of  two 


parts.  One  part,  is  akin  to  the  response  of  the  uncoupled 

rail-fastener  system  to  a point  force  excitation  and  the  other 
part,  Ig,  closely  resembles  the  response  of  the  uncoupled  slab 
system.  In  w,  these  two  parts  make  a nearly  equal  contribution, 
while  C consists  mainly  of  I^. 

is  little  affected  by  damping  (which  serves  only 
to  attenuate  the  propagating  bending  waves)  except  when  o)  is 
less  than  a damping  bandwidth  away  from  o)  = oj^^/1  + ~ ^or’ 

Near  this  frequency,  the  amplitude  of  1^,  rises  to  a sharp  maximum 
limited  only  by  the  damping.  Similar  remarks  apply  to  1^,  except 
that  its  maximum  occurs  at 


Thus  both  Icl  and  Iw,  i show  resonant  maxima  near  o)  and  at  ui 

I ^ ' 1 ' rvr  nc: 


OS 


Another  topic  of  interest  here  is  the  behavior  at 
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structural  wave  coincidence.  Strictly  speaking,  structural 
wave  coincidence  occurs  when  the  dispersion  curve  of  the  rail 
motion  with  slab  displacement  blocked  crosses  the  dispersion 
curve  of  the  slab  motion  with  the  rail  displacement  blocked. 

The  resulting  condition  for  coincidence  is 

k^  = or  A = 0 (C.16) 

For  3 small,  this  means  that  an  approximate  coincidence  condi- 
tion is  simply  the  intersection  of  the  uncoupled  rail  and  slab 
dispersion  curves  determined  above. 

It  is  instructive  to  tie  the  coincidence  condition  to 
the  general  behavior  of  the  coupled  system  mode  dispersion 
curves.  These  are  given  by  the  zeroes  of  V in  Eq.  C.IO  in  the 
absence  of  damping.  Denote  the  coincidence  wavenumber  and  fre- 
quency by  k and  oo  . If  we  consider  wavenumbers  and  frequencies 

far  from  k and  w so  that  |A|  >>  k k gijj,  then  X in  Eq.  C.ll  is 
3 

of  order  k^K2^/^  • Thus,  in  this  case  and  for  small  3,  the  zeroes 

4 4 4 4 

of  Eq.  C.ll  are  k k and  k - a , i.e.,  they  are  on  the  un- 

X 3T  X 

4 4 

coupled  system  dispersion  curves.  The  branch  k^  =;  k^  occurring 

4 4 . ^ • 

for  0)  < 0)  becomes  k - a for  w > co  and  vice-versa.  If  coin- 
c X c 

cidence  occurs  at  all,  A must  achieve  a minimum  for  some  k^  and  w. 
From  Eq.  C.ll  this  occurs  at  A = 0;  precisely  the  point  of  coin- 
cidence. This  condition  also  implies  the  minimum  separation  (along 
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k ) of  the  coupled  system  dispersion  curves,  equal  to 


A - 2A  1 . = 2|X  I.  - = 

r'min  'r'A=0  12^1 


(C.17) 


It  then  follows  that  the  total  kinetic  energy  of  rail 

f ''2  ( ''2 

or  slab  (related  to  the  integrals  |dk^|H^j^|  , Jdk^  | | is 

greatest  at  the  coincidence  frequency,  since  here  the  poles  of 
or|H^|^  are  closest  together. 


Finally,  we  give  the  form  of  5 and  w^  at  coincidence. 

For  the  present  calculation,  we  assume  that  I | and  >> 

and  Ag.  Then  at  coincidence: 


A = *s  = ^2  2 


OS 


(C.18a) 


or 


^r  = - ^1 


(C.18b) 


and,  from  Eq.  C.15 


? = 1 1 
- 4 I 1 3 
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A X 
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4 k 


+ j e 


x(j  |k^|  - ^ A^/|k^ 


cos  X 
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1^1 

4|k_  I 


3 ^rl^rl\' 
j ^ ^ 7 's  (C.19a) 

l^rl'  j\ 
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o > 
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1 1 1 f ■''1^ 


4 K 


1 lk 


! 

r ' i. 


, ■ . I l^ri- 

/e  sinh  ; 


\4|k 


x(j  |kj.|  - l!Aj./|kj.|^)  / |X^| 

-e  sin  X 


3 ^16 


3 r ' r 


|k 


(C.19b) 


r ' / 


to  within  terms  of  order  |X^|/|k^l  . 
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APPENDIX  D : REPORT  OF  INVENTIONS 


This  report  presents  the  development  of  a model 
to  predict  the  performance  of  floating  slab  tracks.  After 
a diligent  review  of  the  work  performed  under  this  contract, 
it  was  found  that  no  new  inventions,  discoveries,  or  improve- 
ments of  inventions  were  made. 
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